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Summary 


During  the  three  years  that  the  Superconducting  and  Magnetically  Ordered  Systems  Research 
Program  at  Virginia  State  University  was  supported  by  the  Air  Force  Office  of  Scientific  Research, 
considerable  progress  was  made  in  these  areas  of  study,  and  the  groundwork  was  laid  for  further 
research  that  has  been  extremely  fruitful. 

During  these  three  years,  members  of  the  group,  primarily  Dr,  David  R.  Noakes,  performed  a 
number  of  studies  on  magnetically-ordered  systems.  These  include  f-electron  magnetic  systems, 
Kondo-lattice  CeTSn  (T  =  transition  metal),  strong-collision  dynamics  in  PrP,^,  rare-earth  - 
transition-metal  intermetallics,  non-Feimi  liquid  Ce7Ni3,  range-correlated  magnetic-moment 
variation,  rare-earth  quasicrystals,  neutron  scattering  studies  of  chromium  alloys,  and  magnetic 
frustration  and  competing  exchange.  This  report  will  focus  on  range-correlated  magnetic-moment 
variation,  rare-earth  quasicrystals,  and  magnetic  frustration. 

During  this  same  period,  members  of  the  group,  primarily  Dr.  Carey  E.  Stronach,  performed  a 
number  of  muon-spin-rotation  experiments  on  high-temperature  superconductors  at  the  TRIUMF 
cyclotron  laboratory.  A  major  focus  of  this  work  was  the  interplay  between  superconductivity  and 
magnetic  ordering  in  these  ceramic  materials.  These  phenomena  are  elaborated  upon  in  this  report. 

This  report  also  contains  a  listing  of  the  26  papers  published  by  the  VSU  group  during  this  three- 
year  period.  Finally,  this  program  provided  excellent  research  experiences  for  participating  students, 
most  of  whom  are  members  of  minority  groups. 
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REPRODUCED  FROM 
BEST  AVAILABLE  COPY 


High-Temperature  Superconductivity 

An  important  discovery  in  this  area  was  the  result  of  a  study  of  the  ceramic  compound 
SrjYRui.yCUyOg,  done  at  TRIUMF  using  the  muon  spin  rotation  (/^SR)  technique.  Two  signals  were 
observed  below  the  Curie  and  Neel  temperatures  (both  ~  30K),  one  that  we  associate  with  the 
YRUO4  layers  that  indicates  static  magnetic  ordering,  the  other  that  we  associate  with  the  SrO  layers, 
which  exhibits  superconducting  behavior  in  transverse  field.  These  results  support  the  hypothesis 
that  the  hole  condensate  is  located  in  the  reservoir  layers,  namely  the  SrO  layers. 

These  results  also  imply  a  picture  of  this  compound,  in  which  there  are  alternating  layers  of 
superconducting  and  magnetically-ordered  material,  with  the  magnetic  field  lying  along  the 
respective  plane  faces,  and  alternating  in  direction.  These  results  are  described  in  detail  in  reference 
16,  by  D.R.  Harshman  et  al.  This  paper  was  presented  by  C.E.  Stronach  at  the  Eighth  International 
Conference  on  Muon  Spin  Rotation/Relaxation/Resonance  in  August/September  1999,  held  in  Les 
Diablerets,  Switzerland.  It  is  included  as  an  appendix  to  this  report.  (D.R.  Harshman  et  al,  Physica 
B  289=290, 360  (2000)). 

In  addition  to  the  new  and  exciting  results  described  above,  several  additional  noteworthy  results 
were  obtained,  submitted  or  published  during  this  reporting  period.  A  paper  on  Coexistence  of 
ferromagnetism  and  superconductivity  in  the  hybrid  ruthenate-cuprate  compound  RuSr2GdCu20g 
by  muon  spin  rotation  and  dc  magnetization  (C.  Bernhard  et  al,  collaboration  with  the  Nffl- 
Stuttgart  and  Konstanz  group,  Phys.  Rev.  B  59  (1999)  14,099)  was  published  in  Jrme  1999  that 
shows  manifestations  of  this  coexistence  in  p'^SR  data.  The  zero-field  measurements  indicate  that 
the  ferromagnetic  phase  is  homogeneous  on  a  microscopic  scale  and  accounts  for  most  of  the  sample 
volume.  Also,  the  magnetic  order  appears  not  to  be  significantly  modified  at  the  onset  of 
superconductivity.  This  paper  is  enclosed  as  an  appendix. 

A  paper  on  Low  temperature  vortex  structures  in  underdoped  Bi2Sr2CaCu20g^g  (T.  Blasius  et  al, 
collaboration  with  the  Konstanz-MPI-Stuttgart  group,  published  in  the  Proceedings  of  the  Eighth 
International  Conference  on  Muon  Spin  Rotation,  Relaxation  &  Resonance,  in  Physica  B)  describes 
measurements  on  the  vortex  state  of  Bi-2212  single  crystals.  We  observed  a  reorganization  of  the 
vortex  structure  with  increasing  vortex  density  that  appears  to  be  related  to  a  disorder-induced 
destruction  of  the  vortex  lattice.  Also,  we  found  evidence  for  a  partial  restoration  of  the  vortex  lines 
at  higher  fields.  This  may  be  related  to  the  increase  of  the  vortex-vortex  interaction  and  the 
subsequent  change  fi'om  single  vortex  pinning  to  collective  pirming. 

A  paper  on  Melting  and  dimensionality  of  the  vortex  lattice  in  underdoped  YBa2CugOg  gQ  was 
published  in  Physical  Review  B.  Stronach  and  Noakes  collaborated  with  J.E.  Sonier  of  Los  Alamos 
National  Laboratory,  and  J.H.  Brewer  and  R.F.  Kiefl  and  associates  of  the  University  of  British 
Columbia  in  this  study,  which  revealed  a  vortex-lattice  melting  transition  at  much  lower  temperature 
than  in  the  fully  oxygenated  material.  This  suggests  that  adjacent  layers  of  "pancake"  vortex 
decouple  in  the  liquid  phase.  In  addition,  evidence  is  found  for  a  pinning-induced  crossover  from 
a  solid  3D  to  a  quasi-2D  vortex  lattice. 

A  paper  entitled  Zero-field  muon-spin-rotation  study  of  hole  antiferromagnetism  in  low-carrier- 
density  Y].fSaJ3a2Cu30g ,  by  Stronach,  Noakes,  X.  Wan  (former  VSU  physics  graduate  student), 
and  memljers  of  the  Konstanz-MPI-Stuttgart  group  was  included  in  the  1 998  annual  report  as  a  pre- 
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print.  It  has  since  been  revised  and  published  (Physica  C  311  (1999)  19).  Based  on  data  taken  at 
TRIUMF,  we  showed  that  even  in  the  low-carrier-density  regime,  the  long-range  antiferromagnetic 
ordering  is  still  strongly  frustrated  by  the  spin-glass  mechanism  through  calcium  doping.  Apparent 
transitions  were  observed  between  the  frozen-hole  antiferromagnetic  state  and  the  spin-glassy  state. 
These  data  made  it  all  the  more  apparent  that  the  same  magnetic  interactions  hold  in  both  low- 
carrier-density  (Y,Ca)-123  and  (La,Sr)2Cu04. 


Range-Correlated  Magnetic  Moment  Variation  (RCMMV) 

In  previously  published  work,  Noakes  and  collaborators  reported  a  static  ZF  muon  spin  relaxation 
function  with  anomalously  shallow,  or  even  completely  absent,  minimum  before  the  recovery  to  1/3 
(of  initial  asymmetry)  in  the  spin-frozen  states  of  icosahedral  Al-Mn-Si  [1]  quasicrystals  and 
Ce(Ni,Cu)Sn  Kondo-semimetal  alloys,[2,3]  and  had  developed  the  phenomenological  “Gaussian- 
broadened  Gaussian”  static  ZF  relaxation  function  that  fits  those  data.  [4] 

Noakes  performed  numerical  simulations  demonstrating  that  “range-correlated  moment-magnitude 
variation”  (RCMMV)  can  produce  arbitrarily-shallow  static  ZF-pSR  functions.  The  depth  of  the 
minimimi  is  an  (inverse)  indicator  of  the  correlation  length,  from  zero  for  the  standard  Gaussian 
Kubo-Toyabe,  to  infinite  correlation  length  for  monotonic  static  relaxation.  The  property  varying 
with  a  correlation  range  is  surprising,  however:  it  is  the  magnitude  of  the  ions'  local  moments. 
Moment  orientation  was  always  completely  random  (as  appropriate  for  frozen  spin  glass-like  states), 
but  their  magnitude  also  varied,  with  a  correlation  length.  A  paper  describing  these  results  has  now 
appeared,  which  is  included  as  an  appendix  to  this  report.  (D.R.  Noakes,  J.  Phys.  Cond.  Matter  11. 
1589  (1999))  [5].  At  the  1999  Int.  Conf.  on  pSR,  independent  confirmation  of  this  work  was 
provided  by  the  Columbia  group,  who  observed  shallow  static  ZF  muon  spin  relaxation  functions 
in  the  "two-leg  spin-ladder  system"  Sr(Cu,.,tZn^)03,  and  performed  independent  numerical 
simulations  arriving  at  a  RCMMV  model  for  the  system.  [6] 

[1]  D.R.  Noakes,  A.  Ismail,  E.J.  Ansaldo,  J.H.  Brewer,  G.M.  Luke,  P.  Mendels  and  S.J.  Poon,  Phys. 

Lett.  A  199, 107(1995). 

[2]  S.J.  Flaschin,  A.  Kratzer,  F.J.  Burghart,  G.M.  Kalvius,  R.  Wappling,  D.R.  Noakes,  R.  Kadono, 

I.  Watanabe,  T.  Takabatake,  K.  Kobayashi,  G.  Nakamoto  and  H.  Fujii,  J.  Phys.  Cond.  Matter 
8, 6967  (1996). 

[3]  G.M.  Kalvius,  S.J.  Flaschin,  T.  Takabatake,  A.  Kratzer,  R.  Wappling,  D.R.  Noakes,  F.J. 

Burghart,  A.  Briickl,  K.  Neumaier,  K.  Andres,  R.  Kadono,  I.  Watanabe,  K.  Kobayashi,  G. 
Nakamoto  and  H.  Fujii,  Physica  B  230-232. 655  (1997). 

[4]  D.R.  Noakes  and  G.M.  Kalvius,  Phys.  Rev.  B  56, 2352  (1997). 

[5]  D.R.  Noakes,  J.  Phys.  Cond.  Matter  IT,  1589  (1999) 

[6]  M.I.  Larkin  et  al ,  Proc.  1 999  Int.  Conf.  pSR,  Physica  B  289-290. 1 53  (2000). 


4 


^SR  in  Rare-Earth  Quasicrystals 

Quasicrystals  are  materials  that  have  long-range  orientational  order  without  translational  periodicity. 
Only  a  small  number  of  quasicrystalline  materials  display  any  magnetism.  Noakes  and  collaborators 
had  previously  published  on  pSR  in  spin-glassy  icosahedral  Al-Mn-Si.[l]  In  Dec.  1996  and  July 
1997,  they  performed  pSR  measurements  of  icosahedral  Gd-Mg-Zn  and  Tb-Mg-Zn,  both  of  which 
show  spin-glass  behavior  in  bulk  measurements,  with  spin-freezing  temperatures  below  lOK. 

Whereas  neutron  scattering  researchers  had  reported  magnetic  "quasi-lattice"  ordering  at  20K  in  Tb- 
Mg-Zn,[3]  our  pSR  spectra  showed  no  evidence  of  any  transition  around  20K,  seeing  instead  spin 
freezing  typical  of  concentrated-moment  spin  glasses  below  lOK  in  both  samples.  Note  that  no  bulk 
measurement  has  indicated  any  transition  near  20K,  so  the  neutron  scattering  report  stands  isolated 
and  imconfirmed.  The  concentrated-spin-glass  behavior  (temperature-dependent  relaxation  rate  to 
far  above  Tg,  zero-field  (ZF)  relaxation  rate  exceeding  instrument  resolution  below  Tg)  for  8  magnetic 
ions  per  hundred  (RE8Mg42Zn5o)  contrasts  with  dilute  spin-glass  behavior  (temperature-dependent 
relaxation  only  below  ~2Tg,  ZF  relaxation  resolved  at  all  temperatures)  in  Aljj  5Mn2o,5Si4  quasicrystal, 
where  there  are  20.5  Mn  ions  per  hundred.  This  re-enforces  the  conclusion  of  the  earlier  work  that 
only  a  fraction  of  the  Mn  ions  are  magnetic,  and  indicates  that  the  magnetism  of  the  rare  earth 
quasicrystals  is  less  anomalous  than  that  of  /- Al-Mn-Si.  Much  faster  relaxation  in  the  Tb  quasicrystal 
Aan  in  the  Gd  one  indicates  significant  "Crystalline  Electric  Field"  (CEF)  effects  for  the  non-S  rare 
earth  ions  in  this  structure,  but  we  did  not  understand  at  that  time  why  power-exponential  (exp[- At]’’, 
where  A  and  p  are  temperature  dependent)  relaxation,  often  foimd  in  pSR  of  concentrated-moment 
spin  glasses  above  Tg,  occurs  only  in  the  Tb  sample  here,  not  the  Gd  one.  This  question  was  resolved 
in  2002. 

While  the  ZF-pSR  became  too  fast  to  measure  below  Tg  in  both  samples,  we  did  observe 
longitudinal-field  (LF)  partial  decoupling  in  applied  fields  up  to  2  Tesla  at  dilution  refrigerator 
temperatures,  and  were  able  to  fit  the  spectra  to  a  model  of  a  very-nearly  static  Gaussian  field 
distribution,  consistent  with  our  dense  spin  glass  description.  These  results  were  published  in  Physics 
Letters  A.  [2]  This  paper  is  included  as  an  appendix  to  this  report. 

[1]  D.R.  Noakes,  A.  Ismail,  E.J.  Ansaldo,  J.H.  Brewer,  G.M.  Luke,  P.  Mendels  and  S.J.  Poon,  Phys. 

Lett.  A 199, 107(1995). 

[2]  D.R  Noakes,  G.M.  Kalvius,  R.  Wappling,  C.E.  Stronach,  M.F.  White,  H.  Saito  and  K. 

Fukamichi,  Phys.  Lett.  A  238, 197  (1998). 

[3]  B.  Charrier,  B.  Ouladdiaf  and  D.  Schmitt,  Phys.  Rev.  Lett.  78, 4637  (1997). 
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Magnetic  Frustration  and  Competing  Exchange 

Frustration  of,  or  competition  between,  magnetic  ordering  interactions  can  reduce  ordering 
temperatures,  can  cause  complicated  magnetic  ordering  structures,  and  can  cause  anomalous 
magnetic  dynamics  above  the  ordering  temperatures.  In  beginning  to  study  these  effects  with  pSR, 
we  have  done  experiments  on: 

GdMm.  where  the  Mn  sub-lattice  is  frustrated  but  the  Gd  lattice  is  not.  Our  pSR  is  providing  new 
information  on  the  several  ordered  states,[l]  which  are  not  well  imderstood,  partly  because  Gd 
absorbs  so  many  neutrons  that  neutron  diffraction  is  difficult  to  do. 

YMni;.  a  transition-moment-only  baseline  material  for  comparison  to  REMn,2,  REFe4Alg  and 
REFcjAlg  competing-exchange  materials,  which  all  have  in  fact  the  same  crystal  structure.  In  spite 
of  three  crystallographically-distinct  Mn  sites  with  two  different  moments  in  the  ordered  state,[2]  our 
paper  presented  at  the  1999  Int.  Conf.  pSR  [3]  shows  that  well  below  ZF  muon  spin  relaxation 
is  fairly  slow  and  characteristic  of  a  Lorentzian  field  distribution  (there  is  no  coherent  oscillation), 
indicating  near-perfect  cancellation  of  all  local  fields  at  the  muon.  Noakes  has  used  his  muon  site- 
search  software  [4]  to  identify  the  likely  interstitial  muon  sites  as  the  xmit  cell  positions  (0,/4,‘A)  and 
(0,0, /4).  These  are  both  high-symmetry  sites  for  which  Noakes  has  now  demonstrated  that  the  dipole 
field  sum  at  both  these  sites  for  the  published  magnetic  structure  is  zero.  With  this  site  identification, 
pSR  in  YMn,2  is  understood.  The  site  identification  extends  to  the  competing-exchange  materials 
with  magnetic  rare  earths,  which  should  make  interpretation  of  pSR  in  them  easier,  as  well. 

TbFegAlfi  and  ErFe^Alg.  c?-band  vs.  RKKY  competing  exchange  materials.  The  very  rapid  relaxation 
caused  by  the  large  rare  earth  moments  limits  the  amoimt  of  information  that  can  be  extracted  reliably 
from  pSR,  but  there  are  clear  indications  of  substantial  spin  dynamics  affecting  the  muon  relaxation 
in  the  wide  range  of  temperature  where  bulk  magnetization  indicates  that  the  ferrimagnetic  moment 
does  not  develop  in  a  standard  Brillouin-curve  manner.  [5] 

TbCo^Ni,  which  is  a  competing  exchange  system  in  the  sense  that  the  magnetic  ordering  switches 

from  antiferromagnetic  for  Ni-rich  to  ferromagnetic  for  Co-rich  but  the  transition  ions  themselves 
are  non-magnetic,  so  there  is  only  the  RKKY  interaction.  For  intermediate  values  of  x  there  are  two 
magnetic  transitions  as  temperature  is  lowered,  with  a  modulated  antiferromagnetic  structure  at  the 
lowest  temperatures.  In  the  intermediate-x  alloys,  pSR  sees  considerable  local  spin  disorder  and 
dynamics  at  all  temperatures  and  in  all  the  magnetic  states.  [6] 

[1]  E.M.  Martin,  E.  Schreier,  G.M.  Kalvius,  A.  Kratzer,  O.  Hartmann,  R.  Wappling,  D.R.  Noakes, 

K.  Krop,  R.  Ballou  and  J.  Deportes,  Physica  B  289-290. 265  (2000). 

[2]  J.  Deportes  and  D.  Givprd,  Solid  State  Commun.  19,  845  (1976). 

[3]  G.M.  Kalvius,  K.H.  Munch,  A.  Kratzer,  R.  WSppling,  D.R.  Noakes,  J.  Deportes  and  R.  Ballou, 

1999  International  Conference  on  Muon  Spin  Relaxation,  Physica  B.289-290. 261  (2000). 

[4]  D.R.  Noakes,  J.H.  Brewer,  D.R.  Harshman,  E.J.  Ansaldo,  and  C.Y.Huang,  Phys.  Rev.  B  35, 6597 

(1987). 

[5]  G.M.  Kalvius,  D.R.  Noakes,  G.  Grosse,  W.  Schafer,  W.  Kockelmann,  S.  Fredo,  I.  Halevy  and  J. 

Gal,  Physica  B  289-290. 225  (2000). 

[6]  G.M.  Kalvius,  D.R.  Noakes,  R.  Wappling,  G.  Grosse,  W.  Schafer,  W.  Kockelmann,  J.K. 
Jakinthos,  and  P.A.  Kotsanides,  Physica  B  326. 465  (2003). 
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Mechanical  Milling 

This  project,  as  originally  proposed  to  the  Air  Force  Office  of  Scientific  Research,  contained  a 
mechanical-milling  component,  which  involved  development  of  a  laboratory  at  Virginia  State 
University  to  produce  and  characterize  nanostructured  iron  alloys.  Dr.  Anthony  S.  Arrott,  Professor 
Emeritus  from  Simon  Fraser  University,  joined  the  VSU  faculty  as  a  part-time  Distinguished  Senior 
Research  Professor  to  provide  scientific  direction  for  this  part  of  the  research  program  in  February 
1998.  Shortly  after  this  program  was  funded,  a  separate  grant  from  the  Ballistic  Missile  Defense 
Organization  funded  the  mechanical  milling  program.  Details  of  the  accomplishments  of  this  program 
are  given  in  its  annual  reports,  which  were  submitted  to  the  Air  Force  Office  of  Scientific  Research 
(F49620-97- 1-0532)  and  the  Ballistic  Missile  Defense  Organization. 

This  program  was  continued  with  support  from  another  grant  fi'om  the  AFOSR  as  the  Iron  Alloys 
Research  Program,  through  which  three  major  discoveries  were  subsequently  made  in  2002-03 :  spin- 
density  waves  in  Fe(Al)  alloys,  magnetic  foam,  and  inverse  melting. 


Magnetism  Laboratory 

A  proposal.  Magnetic  Materials  Laboratory  Development,  submitted  to  the  F  Y  99  DOD  HBCU/MI 
Infrastructure  Support  Program,  was  funded  for  $194,540  plus  $5,000  matching  funds  by  Virginia 
State  University,  effective  September  1 , 1 999.  This  grant  has  enabled  the  VSU  researchers  to  develop 
a  magnetometer  system  of  unprecedented  sensitivity.  Equipment  purchased  includes  a  cryofree 
superconducting  split-coil  magnet,  a  Compumotor  linear  servo  drive,  a  sensitive  digital  voltmeter, 
a  differential  scanning  calorimeter,  a  bipolar  power  supply,  and  a  precision  LCR  meter. 

This  laboratory  is  focused  primarily  on  the  fabrication  and  characterization  of  nanostructured  iron 
alloys,  which  is  conducted  by  Dr.  Anthony  S.  Arrott,  Distinguished  Senior  Research  Professor  of 
Physics,  Stronach,  Noakes  and  students. 

In  addition  to  the  equipment  to  be  purchased  through  this  grant,  Arrott  has  obtained,  as  a  gift  from 
the  Honeywell  Corporation,  a  2.0-Tesla  Walker  Scientific  electromagnet,  which  plays  an  important 
role  in  the  magnetic  measurements  planned  for  this  laboratory.  The  faculty  participating  in  this 
project  include  Stronach,  Arrott,  Noakes,  and  George  W.  Henderson,  Assistant  Professor  of  Physics. 

International  Muon  Spin  Rotation  Conference 

As  is  noted  at  several  places  in  this  report,  Stronach  and  Noakes  attended  the  Eighth  International 
Conference  on  Muon  Spin  Rotation/Relaxation/Resonance,  which  was  held  in  Les  Diablerets, 
Switzerland,  August  30  -  September  3, 1 999.  They  presented  a  number  of  papers,  to  which  references 
are  made  in  this  report.  It  should  also  be  noted  that  Stronach,  who  was  a  member  of  the  International 
Advisory  Committee  on  Muon  Spin  Rotation,  issued  an  invitation  to  the  group  to  hold  the  2002 
conference  in  Virginia.  The  invitation  was  accepted,  and  the  precise  location  was  later  chosen  to  be 
Williamsburg.  Stronach  and  Noakes  were  named  chairman  and  secretary,  respectively. 
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Ceasar  U.  Jackson,  Jr.  completed  his  studies  and  received  the  Master  of  Science  degree  in  physics 
at  the  May  2000  commencement  exercises  at  Virginia  State  University.  His  thesis  is  titled  "Muon 
Spin  Rotation  Study  of  the  Anomalous  /-electron  Compoimd  CcyNij". 

Morris  F.  White,  Jr.  completed  his  thesis,  "Numerical  Investigations  of  Novel  Observed  Muon  Spin 
Relaxation  Functions",  in  December  2000,  and  received  his  Master  of  Science  degree  in  May  2001 . 

Both  Jackson  and  White  visited  the  TRIUMF  laboratory  twice  to  obtain  data  for  their  theses.  White's 
travel  expenses  were  covered  by  his  employer,  Philip  Morris,  USA. 

Three  undergraduates,  Larry  Harris,  Jerrod  Spencer,  and  Mensah  Alkebu-lan  participated  in  the 
research  program  and  graduated  with  Bachelor  of  Science  degrees  in  physics.  Harris  then  attended 
Clemson  University,  where  he  received  his  Master  of  Science  degree  in  May  2003.  Spencer  was 
tragically  killed  in  an  automobile  accident  in  early  2003. 


Appendices 

Four  of  the  26  papers  published  with  support  from  AFOSR  grant  F49620-97- 1-0297  are  attached  as 
appendices  to  this  report.  They  were  chosen  because  they  represent  the  scientific  findings  with  the 
greatest  impact  that  this  group  published  as  a  direct  consequence  of  this  support.  The  term  "direct 
consequence"  is  used  because  this  support  for  the  magnetism  research  team  at  Virginia  State 
University  has  led  to  many  additional  important  findings,  and  we  are  confident  that  many  more  are 
yet  to  come. 

The  attached  papers  are; 

[1]  D.R.  Harshman,  W.J.  Kossler,  A.J.  Greer,  C.E.  Stronach,  E.  Koster,  B.  Hitti,  M.K.  Wu,  D.Y. 
Chen,  F.Z.  Chien,  H.A.  Blackstead,  J.D.  Dow,  Location  of  the  superconducting  hole  condensate  in 
Sr2YRu,.^CujOs  by  ffSR,  Physica  B  289-290,  360  (2000). 

[2]  C.  Bernhard,  J.L.  Tallon,  C.  Niedermayer,  T.  Blasius,  A.  Golnik,  E.  Brucher,  R.K.  Kremer,  D.R. 
Noakes,  C.E.  Stronach  and  E.J.  Ansaldo,  Coexistence  of  ferromagnetism  and  superconductivity  in 
the  hybrid  ruthenate-cuprate  compound  RuSr2GdCu20g  studied  by  muon  spin  rotation  and  DC 
magnetization.  Physical  Review  B  59, 14099  (1999). 

[3]  D.R.  Noakes,  A  correlation  length  measured  by  zero-field  muon  spin  relaxation  in  disordered 
magnets.  Journal  of  Physics:  Condensed  Matter  11, 1589  (1999). 

[41  D.R  Noakes,  G.M.  Kalvius,  R.  Wappling,  C.E.  Stronach,  M.F.  White,  H.  Saito  and  K. 
Fukamichi,  Spin  dynamics  and  freezing  in  magnetic  rare-earth  quasicrystals,  Phys.  Lett.  A  238. 1 97 
(1998). 
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G.M.  Kalvius,  D.R.  Noakes,  and  O.  Hartmann. 
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Abstract 

The  magnetic  and  superconducting  characteristics  of  sintered  Sr2  YRu^  .^Cu^Og  (for  u  =  0.05, 0.10, 0.15)  were  probed 
using  transverse-  and  zero-field  muon-spin  rotation  (p'^SR).  Since  positive  muons  are  attracted  to  oxygen  ions  in  the 
high-Tc  oxides,  Sr2YRui_uCu„06  should  (and  does)  present  two  types  of  p"^  sites,  one  in  the  YRUO4  layers  and  the 
other  in  the  SrO  layers.  The  transverse-  and  zero-field  pSR  spectra  for  all  three  stoichiometries  exhibit  magnetic  ordering 
at  and  below  w  30  K,  with  a  static  local  field  of  3kG.  This  transition  is  marked  by  a  dramatic  increase  in  the 
p"^  spin  relaxation  rate  as  the  temperature  decreases  below  corresponding  to  an  increased  static  disordering  of  the 
magnetic  moments.  Above  no  static  fields  are  observed.  Instead,  the  data  exhibit  a  slow  dynamic  depolarization, 
presumably  due  to  rapid  fluctuation  of  magnetic  moments.  Both  transverse-  and  zero-field  data  also  indicate  a  smaller 
second  component  ( 10%)  that  we  associate  with  the  SrO  layer,  exhibiting  superconducting  behavior  in  transverse  field 
with  an  observed  bulk  30  K.  ©  2000  Elsevier  Science  B. V.  All  rights  reserved. 

Keywords:  Superconductivity;  Antiferromagnetism;  Muon  spin  rotation 


High-Tc  superconductivity  has  (in  the  past)  gen¬ 
erally  been  associated  yiith  the  perovskite  struc¬ 
tures  containing  CUO2  planes.  However,  it  has 
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been  suggested  that  CUO2  planes  are  not  necessary 
for  high-Tc  superconductivity.  Indeed,  some  or¬ 
ganic  superconductors  are  included  in  the  growing 
list  of  high-Tc  superconductors  [1,2].  In  this 
work,  we  describe  a  detailed  study  of  the  magnetic 
order  and  superconductivity  of  the  Cu-doped 
Sr2YRui_„Cu„06  system  {u=0.05,  0.10,  0.15) 
with  so  little  Cu  that  CUO2  planes  do  not  form 
(see  Fig.  1).  In  the  high-Tc  oxides  the  implanted  \i^ 
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Fig.  1.  On  the  right  is  the  crystal  structure  of 
Sr2YRui_uCu„06.  Two  nearly  equivalent  oxygen  sites  in  the 
YRuj  _„CUu04  layer  are  designated  0(1)  and  0(2),  and  the  0(3) 
site  is  in  the  SrO  layer.  On  the  left  is  the  crystal  structure  of 
La2-„Sr„Cu04  for  comparison. 


particles  tend  to  be  trapped  near  the  oxygen  sites 
[3].  In  the  case  of  Sr2YRu06,  there  are  three  in¬ 
equivalent  oxygen  sites:  Two  in  the  YRUO4  layer 
[designated  0(1,2),  because  the  0(1)  and  0(2)  sites 
are  only  slightly  different],  and  one  in  the  SrO  layer 
[designated  0(3)].  There  are  correspondingly  two 
muon  sites  available.  One  is  formed  by  four  (two 
0(1)  and  two  0(2))  oxygens  in  the  YRui_„Cu„04 
layer,  designated  the  |ao(i,2)-site.  The  second  site, 
associated  with  the  SrO  layer  and  designated  the 
lio(3)-site,  is  also  formed  by  four  oxygens,  two  0(3) 
and  an  0(1)  and  an  0(2).  The  volume  available  to 
the  muon  at  each  of  these  sites  is  roughly  commen¬ 
surate  with  the  size  of  muonium.  Thus  we  expect 
to  see  a  two-component  muon-spin  relaxation 
spectrum. 

Polycrystalline  samples  of  Sr2YRui_uCu„06 
were  prepared  by  the  standard  solid-state  reaction 
method  [4].  Structural  characterization  was  car¬ 
ried  out  with  an  energy  dispersive  X-ray  analyzer, 
by  high-resolution  X-ray  diffraction,  and  by  neu¬ 
tron  scattering  —  confirming  that  the  samples  have 
the  Sr2  YRuOe  stoichiometry  and  phase  homogen- 


Temperature  (K) 

Fig.  2.  The  exponential  relaxation  rate  X  versus  temperature  for 
the  Po(i,2)  (squares)  and  poo)  (circles)  sites.  The  triangles  repres¬ 
ent  data  taken  by  first  zero-field  cooling  and  then  applying  the 
500  G  transverse  field.  The  inset  shows  X{T)  for  the  poo)  site  on 
a  linear  scale.  The  lines  connecting  the  data  are  to  guide  the  eye. 


city.  The  upper  limit  of  the  amount  of  impurity 
phase  is  estimated  to  be  less  than  1%.  The  standard 
time-differential  p'^SR  technique  was  employed, 
utilizing  the  M20  beam  line  at  the  TRIUMF  accel¬ 
erator  [5,6].  Our  spectrometer  incorporates  a  veto 
system  that  virtually  eliminates  background  signals 
from  muons  that  miss  the  sample,  which  greatly 
facilitates  the  recovery  and  separation  of  a  minority 
signal.  In  this  work,  we  utilize  the  local-probe  na¬ 
ture  of  the  p'^'SR  technique  to  obtain  microscopic 
evidence  for  the  existence  of  bulk  superconductivity 
in  this  material,  which  on  the  one  hand  is  magnetic 
and  on  the  other  superconducts  although  it  pos¬ 
sesses  neither  CuO  nor  CUO2  planes. 

Figs.  2  and  3  show  the  transverse-field  relaxation 
rate  X  and  precession  frequency  v  versus  temper¬ 
ature  for  w  =  0.15  (m  =  0.05  and  0.10  show  similar 
behavior).  At  temperatures  below  ^  30  K,  both  the 
zero-field  (ZF)  and  transverse-field  (TF)  spectra 
indicate  two  components,  designated  po(i,2)  and 
Po(3)j  reflecting  the  expected  two  distinct  muon 
sites.  The  subscripts  correspond  to  the  oxygen  in 
the  YRUO4  layers  and  the  oxygen  in  the  SrO  layers, 
respectively.  In  ZF,  the  jUo<i,2)  site  (  90%  of  total 

signal  amplitude)  exhibits  fast  relaxation  and 
a  large  ( 3  kG)  local  magnetic  field,  while  the 
Po(3)  site  (  ~  10%  of  total  amplitude)  shows  neither 
significant  relaxation  nor  a  locally  ordered  field. 
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Fig.  3.  The  muon  precession  frequency  versus  temperature  for 
the  po(i.2)  (squares)  and  poo)  (inset  -  circles)  sites.  The  triangles 
represent  data  taken  by  zero-field  cooling  and  then  applying  the 
500  G  transverse  field.  The  lines  connecting  the  data  points  are 
to  guide  the  eye. 


From  these  data,  the  spin-freezing  temperature, 
Tn,  appears  to  be  30  K.  At  low  temperatures, 
the  3kG  local  field  at  the  |io(i,2)  muon  site  is 
necessarily  associated  with  the  ordering  of  the  Ru 
moments.  Moreover,  the  very  fast  relaxation  of  the 
muon  ensemble  at  the  Po(i,2)  site  reflects  a  strong 
disordering  of  the  Ru  moments  as  well.  When  com¬ 
bined  with  recent  polarized  neutron  data  [7],  a  pic¬ 
ture  of  in-plane  freezing  of  Ru  spins  with  the 
average  polarization  reversing  direction  from  one 
YRUO4  layer  to  the  next  emerges,  thereby  presenti¬ 
ng  a  zero  field  at  the  poo)  site  in  the  SrO  layer. 
Above  30  K  relaxation  due  to  static  electronic 
moments  disappears,  suggesting  motional  narrow¬ 
ing  of  the  lines  due  to  rapid  fluctuations  of  the  Ru 
moments. 

Upon  application  of  a  500  G  transverse  field,  the 
Po(i,2)  site  again  displays  rapid  spin  depolariz¬ 
ation  and  large  local  magnetic  field,  while  the 
Po(3)  site  clearly  exhibits  4  10-fold  rise  in  the  relax¬ 
ation  rate  X  as  temperature  decreases  from  30  to 
2K,  along  with  a  concomitant  diamagnetic  shift 
-  consistent  with  superconductivity.  Exponential 
depolarization  forms  [exp(  —  Xt)]  were  assumed  in 
analyzing  both  signals.  An  additional  DC  offset 
was  included  since  the  apphed  field  was  much 
smaller  than  the  local  field.  Note  the  relative  in¬ 


crease  in  A  at  2  K  for  the  poo)  site  after  zero-field 
cooling  (ZFC  -  see  inset  of  Fig.  2).  This  is  an 
expected  consequence  of  an  induced  increase  in  the 
spatial  disorder  of  the  superconducting  vortices. 
Interestingly,  when  the  magnetic  field  is  removed  in 
a  field-cooled  sample,  we  see  little  evidence  of  trap¬ 
ped  flux,  indicating  very  weak  flux  pinning.  This 
may  be  related  to  having  essentially  zero  vortex 
coupling  between  superconducting  layers,  i.e.,  iso¬ 
lated  sheets  of  “pancake”  vortices,  which  is 
consistent  with  having  the  superconducting  hole 
condensate  residing  in  the  SrO  layers.  Thus,  fully 
developed  bulk  superconductivity  forms  below 
Tn  ^  30  K  (i.e.,  when  static  magnetic  order  is  seen 
at  the  Po(i,2)  site). 

We  now  compare  our  data  with  two  competing 
models.  The  dominant  fact  about  this  material  is 
that  superconductivity  and  magnetism  co-exist  in 
the  same  unit  cell  (not  necessarily  in  the  same 
layer).  The  Ren-Wu  p-wave  pairing  model  [8-10] 
has  ferromagnetic  coupling  by  double  exchange 
between  adjacent  Ru  spins  [11,12],  and  was  solved 
in  a  molecular-field  approximation  by  de  Gennes 
[13,14].  The  magnetic  properties  are  determined  by 
J,  the  AFM  exchange,  and  fc,  the  inter-layer  double 
exchange.  Monte  Carlo  calculations  revealed  that, 
for  a  given  strength  of  the  interlayer  double  ex¬ 
change,  two  magnetically  ordered  states  can  be 
encountered  as  the  temperature  is  lowered:  (i)  anti- 
ferromagnetic  order  (for  Ti  <  T  <T^)  followed 
by  canting  (for  T  <  Ti),  or  (ii)  ferromagnetic  order 
(for  Ti  <  T  <  Tp  )  followed  by  canting  (for 
T  <  Ti).  Here,  Ti  is  the  spin-canting  temperature, 
Tp  is  the  ferromagnetic  order  temperature,  and 
Tn  is  the  Neel  temperature.  Significant  predictions 
of  this  picture  are  (i)  that  the  superconductivity  (SC) 
will  have  odd  parity  and  p-wave  symmetry,  and  (ii) 
the  current-carrying  holes  of  SC  will  reside  prim¬ 
arily  in  the  YRUO4  layers.  However,  the  requisite 
canting  behavior  has  not  been  observed  [7];  this 
picture  remains  unsupported  by  experiments. 

The  charge-reservoir  oxygen  model  (see  e.g.  Ref. 
[15-20]),  on  the  other  hand,  assigns  the  hole  con¬ 
densate  to  the  reservoir  layers,  namely  the  SrO 
layers  in  these  materials.  Its  pairing  is  conventional 
BCS-like  (in  all  likelihood,  s-wave  or  extended 
s-wave),  but  allows  for  electronic  as  well  as  phononic 
pairing.  It  also  embraces  the  bond-valence-sum 
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calculations  [21-25]  for  the  ionic  charges,  which 
put  the  holes  in  the  SrO  layers.  Significant  predic¬ 
tions  of  this  model  include:  (i)  provided  that  Gd 
does  not  order  at  a  high  temperature, 
BaaGdRui  _„Cu„06  should  not  superconduct,  due 
to  pair-breaking  by  paramagnetic  L  =  0  Gd;  (ii)  the 
under-charged  oxygen  ions  that  carry  the  primary 
supercurrent  are  in  the  SrO  layers,  where  there  is 
no  large  average  magnetic  field;  (iii)  the  |x''’SR  sig¬ 
nal  amplitude  from  the  SrO  layers  will  be  weaker 
than  the  signal  amplitude  from  the  YRUO4  layers, 
because  the  holes  (repulsive  to  n'*')  reside  (prim¬ 
arily)  in  the  SrO  layers  where  the  oxygen  ions  are 
not  fully  charged  to  O”^  [15-20];  and  (iv)  the 
carriers  of  SC  are  holes  and  so  the  Hall  coefficient 
should  be  positive  (hole-like).  The  present  results 
support  this  model. 

In  summary,  our  data  indicate  that  the 
Po(i,2)  site  should  be  associated  with  oxygens  in 
the  YRUO4  layers  and  the  poo)  site  with  the  SrO 
layers.  The  marked  difference  in  the  trapping 
probabilities  between  the  two  sites  can  be  under¬ 
stood  as  due  to  a  reduced  negative  charge  on  the 
SrO  oxygen:  this  is  where  the  current-carrying 
holes  reside.  Muons  stopped  near  oxygen  ions  in 
the  SrO  layers  give  evidence  for  vortices  below 
~  30K,  as  revealed  by  enhanced  disorder  seen  in 
the  TF  data  after  ZFC  at  low  temperature.  We 
suggest  that  it  may  not  be  a  coincidence  that  the 
in-plane  static  magnetic  ordering  of  the  Ru  mo¬ 
ments  and  the  observed  fully  developed  supercon¬ 
ductivity  in  the  SrO  layers  both  occur  at  ~  30  K. 
We  further  posit  that  the  intrinsic  Tc  (i.e.,  in 
the  absence  of  fluctuating  magnetic  moments) 
of  the  superconducting  condensate  may  be  higher, 
but  the  rapid  fluctuations  of  the  Ru  moments  above 
Tn  may  suppress  the  superconducting  state.  Note 
that  microwave  measurements  reported  elsewhere 
[26]  ‘  indicate  flux  dissipation  above  30  K.  When 
the  temperature  is  decreased  below  T^,  the  fluctu¬ 
ations  cease  and  the  in-plane  magnetic  ordering 
(with  average  polarization  reversing  direction  from 
one  YRUO4  layer  to  the  next)  is  established.  This 
unique  magnetic  state  (which  provides  a  zero  net 


'  Microwave  measurements  on  the  same  samples  indicate  an 
onset  Tc  of  42.5  K. 


magnetic  field  in  the  SrO  layers)  may  no  longer 
inhibit  hole  pairing,  but  does  result  in  the  observed 
fully  developed  superconductivity  in  the  SrO 
layers.  Our  results,  along  with  recent  polarized 
neutron  scattering  data  [7],  are  consistent  with  a 
ferromagnetic  YRui_„Cu„04  layer  (the  data  show 
a  broadened  field  distribution)  in  which  the  mag¬ 
netization  vector  lies  in  the  a-b  plane  and  reverses 
direction  from  magnetic  layer  to  magnetic  layer 
producing  a  net  antiferromagnetic  ordering.  While 
these  results  are  completely  consistent  with  the 
predictions  of  the  charge-reservoir  oxygen  model 
[15-20],  we  cannot  absolutely  rule  out  the  possibil¬ 
ity  of  exotic  superconductivity  with  higher-order 
pairing  in  the  YRUO4  layers. 
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We  have  investigated  the  magnetic  and  the  superconducting  properties  of  the  hybrid  ruthenate-cnprate 
compound  RuSr2GdCu208  by  means  of  zero-field  muon-spin  rotation  (ZF-/iSR)  and  dc  magnetization  mea¬ 
surements.  The  dc-magnetization  data  established  that  this  material  exhibits  ferromagnetic  order  of  the  Ru 
moments  [/^(Ru)*^!  fig]  below  rc=133  K.  and  becomes  superconducting  at  a  much  lower  temperature  Tc 
=  16  K.  The  ZF-/xSR  experiments  indicate  that  the  ferromagnetic  phase  is  homogeneous  on  a  microscopic 
scale  and  accounts  for  most  of  the  sample  volrune.  They  also  suggest  that  the  magnetic  order  is  not  signifi¬ 
cantly  modified  at  the  onset  of  superconductivity.  [S0163-1829(99)07321-X] 


I.  INTRODUCTION 

Since  the  discovery  of  superconductivity  in  the  cuprate 
system  La2-.^Ba,Cu04  in  1986,^  an  ever  growing  variety  of 
high-Tc  superconducting  cuprate  compounds  has  been  syn¬ 
thesized  all  of  which  contain  Cu02  planes  (some  also  contain 
CuO  chains)  as  their  essential  structural  elements  which  host 
the  superconducting  charge  carriers.^  Between  the  Cu02 
planes  are  various  kinds  of  layers,  typically  NaCl-type, 
which  are  insulating  and  act  merely  as  a  charge  reservoir.  To 
date,  the  ruthenate  compound  Sr2Ru04  is  the  only  known 
layered  perovskitelike  system  which  becomes  superconduct¬ 
ing  even  though  it  does  not  contain  any  Cu02  planes  or  CuO 
chains.^  Despite  its  rather  low  transition  temperature  Tc 
=  1.5  K,  the  study  of  its  electronic  and  magnetic  properties 
has  become  a  very  rich  and  active  field  of  research.^  In  par¬ 
allel,  the  electronic  and  magnetic  properties  of  the  related 
ruthenate  compounds,  such  as,  for  example,  the  SrRu03  sys¬ 
tem  which  is  an  itinerant  4df-band  ferromagnet  with  Tc 
« 165  K,  have  attracted  a  great  deal  of  interest.^ 

Another  potentially  promising  and  exciting  direction  of 
research  has  been  prompted  by  the  circumstance  that  the 
Ru02  layers  share  the  same  square-planar  coordination  and  a 
rather  similar  bond  length  with  their  Cu02  counterparts.  A 
whole  new  femily  of  hybrid  ruthenate-cuprate  compounds 
may  therefore  be  constructed  whose  members  consist  of  dif¬ 

PRB  59 


ferent  sequences  of  alternating  Ru02  and  Cu02  layers.  Re¬ 
cently,  one  such  a  hybrid  ruthenate-cuprate  compound,  the 
1212-type  RuSr2GdCu208  system  comprising  Cu02  bilayers 
and  Ru02  monolayers,  has  been  synthesized  as  a  single¬ 
phase  material.^  A  subsequent  study  of  its  electronic  and 
magnetic  properties  has  revealed  that  this  material  exhibits 
electronic  ferromagnetic  order  at  a  rather  high  Curie  tem¬ 
perature  Tc- 133- 136  K  and  becomes  superconducting  at  a 
significantly  lower  critical  temperature  7*^=  15-40  K  (de¬ 
pending  on  the  condition  of  preparation  and  annealing). 
The  most  surprising  observation,  however,  is  that  the  ferro¬ 
magnetic  order  does  not  vanish  when  superconductivity  sets 
in  at  .  Instead,  it  appears  that  the  ferromagnetic  state  re¬ 
mains  largely  unchanged  and  coexists  with  superconductiv¬ 
ity.  This  finding  implies  that  the  interaction  between  the  su¬ 
perconducting  and  the  ferromagnetic  order  parameters  is 
very  weak  and  it  raises  the  question  of  whether  both  order 
parameters  coexist  on  a  truly  microscopic  scale.  Since  the 
early  investigations  of  Ginzburg  in  1957,’  the  prevailing 
view  is  that  the  coexistence  of  a  superconducting-  (with  sin¬ 
glet  Cooper  pairs)  and  a  ferromagnetic  order  parameter  is  not 
possible  on  a  microscopic  scale  since  the  electromagnetic 
interaction  and  the  exchange  coupling  lift  the  degeneracy  of 
the  spin-up  and  the  spin-down  partners  of  the  Cooper  pair 
and  cause  strong  pair  breaking.  Indeed,  merely  based  on 
magnetization  and  transport  measurements  one  cannot  ex¬ 
clude  the  possibility  that  the  RuSr2GdCu208  samples  may  be 
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spatially  inhomogeneous  with  some  domains  exhibiting  fer¬ 
romagnetic  order  and  others  superconducting  order.^  We 
note  that  unambiguous  evidence  for  the  occurrence  of  bulk 
superconductivity  in  RuSr2GdCu208  has  recently  been  ob¬ 
tained  from  specific-heat  measurements  which  reveal  a  size¬ 
able  jump  at  Tc  of  A  Cp  IT 0.35  mJ/g  at  character¬ 
istic  of  a  strongly  underdoped  cuprate  superconductor.^®  In 
the  following  we  report  on  muon-spin  rotation  (/^SR)  mea¬ 
surements  which  establish  that  the  ferromagnetic  order  is 
uniform  and  homogeneous  even  on  a  microscopic  scale. 

The  ^tSR  technique  is  ideally  suited  for  such  a  purpose 
since  it  provides  an  extremely  sensitive  local  magnetic  probe 
and,  fiirthermore,  allows  one  to  reliably  obtain  the  volume 
fraction  of  the  magnetically  ordered  phase.  Here  we 
present  the  result  of  a  zero-field  muon-spin  rotation  (ZF- 
)LtSR)  study  of  a  RuSr2GdCu208  sample  with  r^  =  16  K  and 
Tc=  133  K  which  provides  evidence  that  the  magnetic  order 
parameter  is  spatially  homogeneous  and  accounts  for  most  of 
the  sample  volume.  Furthermore,  the  ZF-/zSR  data  establish 
that  the  ferromagnetic  order  is  hardly  affected  by  the  onset  of 
superconductivity  and  persists  to  the  lowest  available  tem¬ 
perature  of  the  experiment  r=2.2  BL  The  ZF-/iSR  data  can 
be  complemented  by  dc-magnetization  measurements  which 
establish  the  presence  of  ferromagnetic  order  from  the  obser¬ 
vation  of  a  spontaneous  magnetization  at  rc=133  K  and  of 
hysteretic  isothermal  magnetic  behavior  with  a  remanent 
magnetization.  It  is  shown  that  the  ferromagnetic  ordering 
involves  the  Ru  magnetic  moments  with  /i(Ru)*«1.05(5)/i5 , 
while  the  larger  Gd  moments  with  ^  (Gd^'^)«^7.4(l)/X5  re¬ 
main  paramagnetic  down  to  very  low  temperatures.  In  addi¬ 
tion,  the  magnetization  measurements  indicate  an  almost 
complete  diamagnetic  shielding  effect  below  Tc . 


IL  EXPERIMENT 

A.  Sample  preparation  and  characterization 

Polycrystalline  samples  of  the  1212-type  system 
RuSr2GaCu20g  have  been  synthesized  as  previously 
described*  by  solid-state  reaction  of  RUO2,  SrC03 ,  Gd203, 
and  CuO  powders.  The  mixture  was  first  decomposed  at 
960  °C  in  air.  It  was  then  ground,  milled,  and  die-pressed 
into  pellets.  The  first  sintering  step  took  place  in  flowing 
nitrogen  atmosphere  at  1010  °C.  This  step  results  in  the  for¬ 
mation  of  a  mixture  of  the  precursor  material  Sr2GdRu06 
and  CU2O  and  is  directed  towards  minimizing  the  formation 
of  SrRu03.^  The  material  was  then  reground  before  it  was 
reacted  in  flowing  oxygen  for  10  h  at  1050  °C.  This  sintering 
step  was  repeated  twice  with  intermediate  grinding  and  mill¬ 
ing.  Each  reaction  step  was  carried  out  on  a  MgO  single¬ 
crystal  substrate  to  prevent  reaction  with  the  alumina  cru¬ 
cible.  Finally  the  samples  were  cooled  slowly  to  room 
temperature  in  flowing  oxygen.  Following  this  procedure  we 
have  also  made  a  Zn-substituted  RuSr2GdCui.942^.06G8 
sample  and  a  Y<-^Gd  cosubstituted  sample 
RuSr2Gdo.9Yo.iCu208.  X-ray-diffiaction  (XRD)  measure¬ 
ments  indicate  that  all  samples  are  single  phase  1212-type 
material  and  give  no  indication  for  traces  of  the  ferromag¬ 
netic  phase  SrRu03.  Figure  1(b)  displays  a  representative 
XRD  spectrum  of  RuSr2GdCu208,  the  plus  signs  show  the 


FIG.  1.  (a)  The  structure  of  RuSr2GdCu208  with  the  Cu  atoms 
sited  at  the  center  of  the  base  of  the  square  pyramids  and  the  Ru 
atoms  at  the  center  of  the  octahedra.  (b)  The  x-ray-diffraction 
(XRD)  spectrum  for  a  RuSr2GdCu208  sample  (Co  K  a  source).  The 
plus  signs  (+)  are  the  raw  x-ray  data  and  the  solid  line  is  the 
calculated  Rietveld  refinement  profile  for  tetragonal  (space  group 
P4/mmm)  RuSr2GdCu2G8. 

raw  data  and  the  solid  line  shows  the  result  of  the  Rietveld 
refinement.  The  related  structure  of  RuSr2GdCu208  is  shown 
in  Fig.  1(a). 

The  electronic  properties  of  RuSr2GdCu208  have  been 
characterized  by  measurements  of  the  temperature-dependent 
resistivity  and  thermoelectric  power.  Representative  results 
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FIG-  2.  (a)  Temperature  dependence  of  the  resistivity  p  of 
RuSr2GdCu20g.  (b)  The  temperature-dependent  thermoelectric 
power  S{T), 

are  shown  in  Figs.  2(a)  and  2(b),  respectively  (see  also  Ref. 
8).  The  temperature  dependence  of  the  thermoelectric  power 
S{T)  and,  in  particular,  its  normal-state  value  of  5(300 
K)*=*75  /xV/K  is  rather  typical  for  a  strongly  underdoped  cu¬ 
prate  superconductor  with  consistent  with  a  hole 

content  of  /7«=*0.07  holes  per  Cu02  planes  and  a  value  of 
^c,max  of  tiie  order  of  100  The  resistivity  measurements 

inicate  that  the  RuSr2GdCu208  sample  exhibits  zero  resis¬ 
tivity  at  a  critical  temperature  of  Tc— 16  K.  The  precise  value 
of  Tc  varies  between  12  and  24  K,  depending  on  synthesis 
conditions,  and  may  be  raised  to  40  K  by  long-term  anneal¬ 
ing.  The  temperature  dependence  of  the  normal-state  resis¬ 
tivity  is  again  characteristic  of  a  strongly  underdoped  super¬ 
conducting  cuprate  compound.  The  ferromagnetic  transition 
at  Tq-  133  K  causes  only  a  small  yet  noticeable  drop  in  the 
resistivity  indicating  that  the  Ru02  layer  is  almost  insulating 
above  Tq  ,  while  being  poorly  conducting  in  the  ferromag¬ 
netic  state.^ 


TIKI 


FIG.  3.  (a)  Temperature  dependence  of  the  zero-field-cooled  dc 
volume  magnetization  RuSr2GdCu2G8.  The  arrows  show  the 
superconducting  and  the  ferromagnetic  transition  at  Tc=\6  IC  and 
Tc=  133  K,  respectively,  (b)  The  field-cooled  molar  magnetization 
Mmol  for  applied  fields  of  H—S.S,  10,  100  Oe. 

compounds,  forms  a  muoxyl  bond  with  one  of  the  oxygen 
atoms.  The  whole  ensemble  of  muons  is  randomly  distrib¬ 
uted  throughout  a  layer  of  100-200  /xm  thickness  and  there¬ 
fore  probes  a  representative  part  of  the  sample  volume.  Each 
muon  spin  processes  in  its  local  magnetic  field  with  a 
precession  frequency  of,  v^=(y;i/2'7r)-B^,  where 
=  135.5  MHz/T  is  the  gyromagnetic  ratio  of  the  positive 
muon.  The  muon  decays  with  a  mean  lifetime  of  r^+ 
«*2.2  /xs“^  into  two  neutrinos  and  a  positron  which  is  pref¬ 
erentially  emitted  along  the  direction  of  the  muon  spin  at  the 
instant  of  decay.  The  time  evolution  of  the  spin  polarization 
P(r)of  the  muon  ensemble  can  therefore  be  obtained  via  the 
time-resolved  detection  of  the  spatial  asymmetry  of  the  de¬ 
cay  positron  emission  rate.  More  details  regarding  the  zero- 
field  (ZF)  /xSR  technique  are  given  below. 


B.  The  technique  of  muon-spin  rotation 

The  muon-spin  rotation  (/^R)  experiments  have  been 
performed  at  the  Ml  5  beamline  of  TRIUMF  in  Vancouver, 
Canada,  which  provides  100%  spin-polarized  muons.  The 
/iSR  technique  is  especially  suited  for  the  study  of  magnetic 
materials  and  allows  one  to  study  the  homogeneity  of  the 
magnetic  state  on  a  microscopic  scale  and  also  to  access  its 
volume  fi:action.^‘  The  /xSR  technique  typically  covers  a 
time  window  of  10“^- 10“’  s  and  allows  one  to  detect  in¬ 
ternal  magnetic  fields  over  a  wide  range  of  0.1  G  to  several 
Tesla.  The  100%  spin-polarized  surface  muons  (E^*«4.2 
MeV)  are  implanted  into  the  bulk  of  the  sample  where  they 
thermalize  very  rapidly  (^^10“^^  s)  without  any  noticeable 
loss  in  their  initial  spin  polarization.  Each  muon  stops  at  a 
well-defined  interstitial  lattice  site  and,  for  the  perovskite 


in.  EXPERIMENTAL  RESULTS 
A.  dc  magnetization 

Before  we  discuss  the  result  of  the  /iSR  experiments,  we 
first  present  some  dc-magnetization  data  which  establish  that 
the  RuSr2GdCu208  sample  exhibits  a  spontaneous  magneti¬ 
zation  at  a  ferromagnetic  transition  of  7^=  133  K  and  be¬ 
comes  superconducting  at  a  much  lower  temperature  T^. 
=  16  K.  Figure  3(a)  shows  the  temperature  dependence  of 
the  volume  susceptibility  xv>  which  has  been  obtained  after 
zero-field  cooling  the  sample  to  7=2  K,  then  applying  an 
external  field  ofi7***=5.5  Oe,  and  subsequently  warming  up 
to  T-  200  K.  The  density  of  the  sample  has  been  assumed  to 
be  p==6.7  g/cm^  corresponding  to  stoichiometric 
RuSr2GdCu208  with  lattice  parameters  of  a =3.84  A  and 
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FIG.  4.  The  isothennal  magnetization  loops  of  RuSr2GdCu20s 
at  r=5,  50,  and  200  K.  The  inset  shows  a  noagnification  of  the 
low-field  region. 

c=  11.57  A.*  The  superconducting  transition  is  evident  in 
Fig.  3(a)  from  the  onset  of  a  pronounced  diamagnetic  shift 
below  7’,,=  16  K,  The  diamagnetic  shift  at  the  lowest  avail¬ 
able  temperature  of  r=  2  K  corresponds  to  an  almost  com¬ 
plete  diamagnetic  shielding  of  the  sample  volume,  implying 
that  at  least  flie  surface  region  of  the  sample  is  homoge¬ 
neously  superconducting.  In  fact,  all  pieces  that  have  been 
cut  from  the  pellet  exhibit  a  similarly  large  diamagnetic 
shielding  effect  (small  differences  can  be  attributed  to  differ¬ 
ent  demagnetization  factors).  Nevertheless,  the  de¬ 
magnetization  measurements  cannot  give  unambiguous  evi¬ 
dence  for  the  presence  of  bulk  superconductivity  since  an 
almost  complete  diamagnetic  shielding  may  also  be  called 
by  a  filamentary  structure  of  superconducting  material  in  a 
small  fraction  of  the  otherwise  nonsuperconducting  material. 
We  note  however,  that  unequivocal  evidence  for  the  occur¬ 
rence  of  bulk  superconductivity  in  RuSr2GdCu20g  has  re¬ 
cently  been  obtained  from  specific-heat  measurements  which 
reveal  a  sizeable  jump  of  A 7=  Cj,/?’®' 0.3 5  irJ/gatK^  at 
Tc,  comparable  to  or  greater  than  that  seen  in  other  under¬ 
doped  cuprates.*®  For  comparison  in  strongly  underdoped 
YBa2Cu307_f  it  is  found  that  A7—O.2-O.3  mJ/gatK^.*"* 
We  also  note  that  the  specific-heat  measurements  have  been 
performed  on  the  same  samples  which  have  been  studied  by 
/iSR-  and  dc-magnetization  measurements.  Figure  3(b)  dis¬ 
plays  the  (low)  field-cooled  molar  magnetization  M„  for  ap¬ 
plied  fields  of  ff“‘=5.5,  50,  and  100  Oe.  The  ferromagnetic 
tr^pgitinn  at  rc=  133  K  is  evident  fiom  the  sudden  onset  of 
a  spontaneous  magnetization.  Evidently,  the  magnetic  order 
parameter  has  at  least  a  sizeable  ferromagnetic  component 
and  it  persists  almost  unchanged  to  the  lowest  measured  tem¬ 
perature  r=2  K.  In  particular,  it  docs  not  appear  to  weaken 
as  superconductivity  sets  in  at  Tc==  16BC.  Additional  evidence 
for  the  presence  of  ferromagnetic  order  is  presented  in  Fig.  4, 
which  shows  that  the  isothermal  magnetization  loops  at  T 
=5  and  50  K  exhibit  hysteretic  magnetic  behavior  with  a 
remanent  magnetization  Oe  at  5  K  and  200  Oe  at 

50  K. 

Having  established  ttie  existence  of  ferromagnetic  order, 
the  question  arises  of  whether  it  involves  the  Ru  moments  or 


FIG.  5.  The  temperature-dependent  inverse  molar  susceptibility 
\lXni  for  the  high-temperature  range  of  400  K>7>100  K.  The  plus 
signs  show  the  best  fit  using  a  two  component  Curie-Weiss  +  Curie 
function.  Shown  in  the  inset  is  the  saturation  magnetization  in  units 
of  effective  Bohr  magnetons  per  imit  volume  as  a  function  of  ap¬ 
plied  field  at  temperatures  of  r=2,  30,  50,  100,  200  K. 

the  Gd  moments.  In  the  following  we  present  high- 
temperature  susceptibility  data  which  indicate  that  the  ferro¬ 
magnetic  order  involves  only  the  Ru  moments,  whereas  the 
Gd  moments  remain  in  the  paramagnetic  state  below  Tc* 
Figure  5  shows  the  inverse  molar  susceptibility,  Mxm 
obtained  for  different  external  fields  in  the 
range  1000  Oe  (solid  lines)  in  the  temperature 

region  200  K<r<400  K.  Shown  by  the  plus  signs  (+)  is  the 
best  fit  to  the  experimental  data  using  a  two-component 
Curie-Weiss + Curie  function,  ;^=Ci/(r-©)  +  C2/7’>  with 
©  =  r^=  133  K  kept  fixed.  This  function  describes  the  ex¬ 
perimental  data  rather  well  and  it  gives  us  very  reasonable 
values  for  the  magnetic  moments,  with  /ii  =  1.05(5  for 
the  moments  that  order  at  and 
moments  that  remain  paramagnetic  below  Tc*.  The  magnetic 
moment  of  the  paramagnetic  component  agrees  reasonably 
well  with  the  expected  magnetic  moment  of  Gd^'^ 
which  for  a  free  Gd^**"  ion^^  is  —  ^nd 

for  the  structurally  similar  GdBa2Cu307-5 
compound.  On  the  other  hand,  the  value  of  the  Ru  mo¬ 
ments  with  /i(Ru)  =  1.05(5)/X5  also  appears  to  be  reason¬ 
able.  For  Ru^*^  the  number  of  Ad  electrons  is  3  and  the 
free-ion  value  of  the  magnetic  moment  is  3  for  the  high- 
spin  state  and  1  for  the  low-spin  state.  The  experimentally 
observed  value  of  /x(Ru)  =  1.05(5  therefore  seems  to 
imply  that  Ru^"^  is  in  the  low-spin  state.  Shown  in  the  inset 
of  Fig.  5  is  the  field-dependent  magnetization  for  the  tem¬ 
peratures  2, 30, 50,  100,  and  300  K.  The  low-temperature 
magnetization  can  be  seen  to  saturate  at  a  value  of 
«8/X5  ,  as  may  be  expected  for  a  system  that  contains  one 
Gd  moment  per  formula  unit  with  /t(Gd)  -1  fiB  P^us  one  Ru 
moment  with  /x(Ru)  =  1  . 

The  idea  that  the  Gd  moments  do  not  participate  in  the 
ferromagnetic  transition  at  7'c=133  K  is  supported  by  the 
result  of  dc-magnetization  measurements  on  the  10%  Y<-+Gd 
cosubstituted  RuSr2Gdo,9Yo,iCu208.  Figures  6(a)  and  6(b) 
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FIG.  7,  The  time-resolved  normalized  mnon-spin  polarization, 
P{t)/P{t=0)y  at  temperatures  of  (a)  r=5.3  K<Tc  and  (b) 
<r=48  K<rc=133  K,  The  large  oscillatory  component  gives 
clear  evidence  for  the  presence  of  a  magnetically  ordered  state. 


FIG.  6.  (a)  The  temperature-dependent  volume  susceptibility  xv 
of  RuSr2GdCu208  (solid  line),  for  Zn-substituted 
RuSr2GdCU|  94Zno.o608  (dotted  line)  and  for  Y<->Gd  cosubstituted 
RuSr2Gdo,9Yo.!Cu208  (dashed  line),  (b)  The  temperature-dependent 
molar  magnetization  shown  by  the  same  symbols  as  in  (a). 

display  the  zero-field-cooled,-  and  the  field-cooled  suscepti¬ 
bilities  (dashed  lines)  and  compare  them  with  the  corre¬ 
sponding  data  on  the  pure  RuSr2GdCu20}{  sample  (solid 
l^e).  It  is  evident  that  the  ferromagnetic  transition  is  not 
significantly  affected  by  the  partial  substitution  of  nonmag¬ 
netic  Y^"**  for  magnetic  Gd^'*'.  Also  shown  in  Figs.  6(a)  and 
6(b)  by  the  dotted  lines  are  the  results  for  the  Zn-substituted 
RuSr2Cui  94Zno.o608  sample.  The  circumstance  that  the  fer¬ 
romagnetic  order  is  not  sdTected  by  the  Zn  substitution  sup¬ 
ports  our  view  that  the  majority  of  the  Zn  impurities  has  been 
introduced  into  the  Cu02  layers  while  hardy  any  of  them 
reside  within  the  Ru02  layers.  Moreover,  we  infer  from  the 
rapid  Tc  suppression  upon  Zn  substitution  that  only  the  Cu02 
layers  host  the  superconducting  charge  carrier  in 
RuSr2GdCu208. 

B.  Zero-field  muon-spin-rotation  (ZF-jLtSR) 

Next  we  discuss  the  result  of  the  zero-field  (ZF)  /iSR 
experiments.  Figures  7(a)  and  7(b)  show  representative  ZF- 
/iSR  spectra  for  the  evolution  of  the  normalized  time- 
resolved  muon-spin  polarization  P(t)IP{0)  at  temperatures 
of  r=5  and  48  K.  The  value  of  the  initial  muon-spin  polar¬ 
ization  P(0)  has  been  determined  by  a  transverse  field  (TF) 
/iSR  experiment  performed  on  the  same  sample  at  a  tem¬ 
perature  above  Tc*  In  the  ferromagnetic  state  below  Tc 
- 133  K  we  find  that  the  spectra  are  well  described  by  the 
relaxation  function: 

P(0/P(^=0)  I  exp(  -Xf)cos(27r( v^)t) 

(1) 


where  (  p^}  is  the  average  muon-spin  precession  frequency 
which  corresponds  to  the  average  value  of  the  spontaneous 
internal  magnetic  field  at  the  muons  sites,  ( 

=  y^/27r(5^),  with  y^=  835.4  MHz/T  being  the  gyromag- 
netic  ratio.  The  damping  rate  of  the  nonoscillating  (longitu¬ 
dinal)  component  A  is  proportional  to  the  dynamic  spin- 
lattice  relaxation  rate  A'^  l/Ti,  whereas  the  relaxation  rate  of 
the  oscillating  (transverse)  component,  X  is  dominated  by  the 
static  distribution  of  the  local  magnetic  field,  i.e.,  X 
Figure  8  shows  the  temperature  dependence  of 
(a)  tile  precession  frequency  (v^)(r),  (b)  the  transverse  re¬ 
laxation  rate  X(r),  and  (c)  the  longitudinal  relaxation  rate 
A{T). 

Before  we  discuss  the  ZF-/iSR  data  in  more  detail,  we 
first  emphasize  the  most  important  implications,  which  are 
evident  from  Figs.  7  and  8.  Firstly,  the  presence  of  an  oscil¬ 
lating  component  in  the  ZF-/iSR  spectra  for  T<Tc— 133  K 
gives  unambiguous  evidence  for  an  ordered  magnetic  state 
which  is  homogeneous  on  a  microscope  length  scale  (of  typi¬ 
cally  20  A).  Secondly,  from  the  amplitude  of  the  oscillating 
component  we  can  deduce  that  the  magnetically 

ordered  state  accounts  for  more  or  less  the  entire  volume  of 
the  sample.  And  thirdly,  from  the  temperature  dependence  of 
the  /xSR  signal  it  becomes  clear  that  the  magnetic  order  per¬ 
sists  almost  unchanged  in  the  superconducting  state. 

L  The  volume  fraction  of  the  magnetic  phase 

In  the  following  we  outline  how  the  volume  fraction  of 
the  magnetically  ordered  phase  is  obtained  from  the  ampli¬ 
tude  of  the  oscillating  component  of  the  ZF-/aSR  spectra.  For 
a  polycrystalline  sample  with  randomly  orientated  grains  in 
zero  external  field  the  local  magnetic  field,  on  average,  is 
parallel  (perpendicular)  to  the  direction  of  the  muon-spin  di¬ 
rection  with  probability  1/3  (2/3).  For  a  homogeneous  mag¬ 
netically  ordered  sample,  one  therefore  expects  that  2/3  of 
the  amplitude  of  the  ZF-/xSR  signal  (the  transverse  compo- 
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FIG.  8.  The  temperature  dependence  of  the  /xSR  signal  of 
RuSr2GdCu208.  (a)  The  muon-spin  precession  frequency,  v^{T) 
(MHz) =135.5  (MHz/T)  (5^),  Shown  by  the  dashed  line  is  the  best 
fit  using  the  scaling  function  v^{T)=vq{\'-TITc^^j  with 
)S=0.333(5),  rc=  133(1)  K,  and  9.70(5)  MHz 

[5^=720(1)G].  (b)  The  relaxation  rate  of  the  precessing  compo¬ 
nent,  X(7^.  (c)  The  relaxation  rate  of  the  nonprecessing  component, 
A(7^'-l/r,, 

nent)  exhibit  an  oscillatory  behavior,  while  1/3  of  the  signal 
(the  longitudinal  component)  is  nonoscillating  and  is  only 
slowly  damped  due  to  spin-flip  excitations.  On  the  other 
hand,  for  a  sample  with  iiihomogeneous  magnetic  order,  for 
example  containing  nonmagnetic  regions,  the  amplitude  of 
the  osciEating  signal  will  be  accordingly  reduced  and  a  sec¬ 
ond  nonoscillating  transverse  component  will  appear.  If  the 
nomhagnetic  regions  are  microscopically  small,  this 
nonosciUating  component  is  Likely  to  have  a  rather  large 
damping  rate  of  the  order  of  due  to  stray  fields 

which  are  imposed  by  the  neighboring  magnetic  domains. 
From  Figs.  7(a)  and  7(b)  it  can  be  seen  that  the  ZF-/iSR  data 
on  RuSr2GdCu208  give  no  indication  for  such  an  inhomoge¬ 
neous  magnetic  state.  As  was  mentioned  above,  the  ampli¬ 
tude  of  the  initial  muon-spin  polarization  P(r=0)  has  been 
determined  from  a  transverse-field  (TF)-/4SR  measurement. 
From  the  size  of  the  amplitude  of  the  oscillatory  component 
we  deduce  that  more  than  80%  of  the  sample  is  magnetically 
ordered  below  Tc=  133  K.  Based  on  this  analysis  we  esti¬ 
mate  that  the  volume  firaction  of  any  disordered  magnetic  or 
nonmagnetic  phase  must  be  well  below  20%.  Note  that  some 
of  the  muons  (typically  10-20%)  do  not  stop  inside  the 
sample  but  somewhere  in  the  cry  estate  walls.  In  the  ZF-/iSR 
experiment  these  muons  give  rise  to  a  missing  fraction  since 
their  spin-polarization  is  much  more  slowly  damped  than  for 


the  rest  of  the  signal.  In  the  TF-/4SR  experiment,  however, 
this  very  slowly  damped  component  can  be  detected  via  its 
precession  in  the  external  field  and  it  contributes  to  the  total 
muon-spin  polarization  P(0).  The  80%  firaction  of  the  mag¬ 
netically  ordered  phase  therefore  has  to  be  regarded  as  a 
lower  bound.  In  fact,  it  is  rather  likely  that  the  entire  sample 
volume  is  magnetically  ordered.  Finally,  we  note  that  the 
muons  apparently  occupy  only  one  muon  site,  since  only  one 
precession  frequency  is  seen  in  the  ZF-/xSR  spectra.  Also  it 
is  clear  from  ^e  ZF-/xSR  data  that  muon  diffusion  effects 
are  negligibly  small  below  rc=133  K,  similar  to  the  other 
cuprate  superconductors  where  muon  diffusion  is  observed 
only  at  significantly  higher  temperatures  of  7>250 

Z  Local  magnetic  field  at  the  muon  site 

It  is  evident  from  Fig.  8(a)  that  the  muon-spin  precession 
frequency  (the  local  field  at  the  muon  site)  does  not  exhibit 
any  strong  anomaly  at  the  superconducting  transition  tem¬ 
perature  Tc .  Instead,  as  shown  by  the  dashed  line,  the  tem¬ 
perature  dependence  of  the  muon-spin  precession  frequency 
(and  thus  of  the  magnetic  order  parameter)  is  well 
described  by  the  function  Vf^(T)  —  VQ(l  —  T/Tc)^y  with  Vq 
=  9.7(1)  MHz  [corresponding  to  (J5^)(r--^0)-720(10)G], 
r^=  133(1)  K,  and  )8=0.333(5).  This  functional  form  is 
strictly  valid  only  in  the  critical  regime  close  to  Tc  but  it  can 
be  seen  to  provide  a  reasonable  description  of  the  magnetic 
order  parameter  over  a  fairly  wide  temperature  range  of  Tc 
^T^  5  K.  The  anomaly  at  very  low  temperature  arises  most 
likely  from  the  magnetic  ordering  transition  of  the  Gd  mo¬ 
ments  at  r^«=*2,6  K.  Note  that  for  the  structurally  related 
compound  GdBa2Cu307_^  (Gd-123)  the  antiferromagnetic 
ordering  transition  of  the  Gd  moment  occurs  at  a  very  similar 
temperature  of  7’i/=  2.3  The  value  of  the  critical  ex¬ 

ponent  /3=0333  is  close  to  the  theoretical  value  0.345  in  the 
3D  XY  model. We  cannot  determine  with  certainty  the 
number  of  components  in  the  spin  system  with  these  data 
and,  in  particular,  distinguish  between  the  two-component 
XY  ()3=0.345)  and  the  three-component  Heisenberg 
=0.365)  models.  The  contribution  of  ferromagnetic  fluctua¬ 
tions  above  Tc  to  the  susceptibility  provides  better  discrimi¬ 
nation  as  will  be  discussed  later. 

The  oscillating  transverse  component  exhibits  a  damping 
rate  of  the  order  of  X*=*10-15  /ts”’  ^corresponding  to  a  spread 
in  the  local  magnetic  field  of  (A5^)/(S^)***^0.2.  This  20% 
spread  of  the  local  magnetic  field  does  not  seem  to  agree 
with  a  scenario  where  the  ferromagnetic  order  is  assumed  to 
exhibit  a  spiral  modulation  (with  a  wavelength  shorter  than 
the  superconducting  coherence  length  of  typically  20  A  in 
the  cuprates)  and/or  to  be  spatially  inhomogeneous  as  in 
ErRh4B4,^’  HoMOfiSg,^®  and  Instead,  we  emphasize 

that  the  observed  spread  in  the  local  magnetic  field  can  be 
accounted  for  by  the  grain-boundary  effects  and  by  the  dif¬ 
ferences  in  the  demagnetization  factors  of  the  individual 
grains  which  naturally  arise  for  a  polycrystalline  sample  that 
has  a  very  small  average  grain  size  of  about  1  /xm.*  Also,  we 
point  out  that  recent  transmission-electron-microscopy  stud¬ 
ies  have  revealed  that  our  present  Ru-1212  sample  contains 
[100]  rotation  tvrins  and  also  exhibits  some  cationic  disorder 
due  to  the  intermixing  of  Sr<->Gd  and  to  a  lesser  extent  of 
Ru4--i^Cu.*  These  kinds  of  structural  imperfections  certainly 
tend  to  further  increase  the  transverse  relaxation  rate  X  of  the 
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ZF-/iSR  spectra.  Meanwhile,  we  have  prepared  Ru-1212 
sanvples  which  are  structurally  more  perfect  (by  sintering  at 
slightly  higher  temperature  and  for  longer  periods).*  Recent 
dc-magnetization  measurements  have  shown  that  these  crys¬ 
tallographic  defects  do  not  affect  the  fundamental  magnetic 
and  superconducting  behavior.  In  &ct,  both  the  supercon¬ 
ducting  and  the  ferromagnetic  transitions  become  somewhat 
sharper  and  T*  and  Tc  are  slightly  increased  for  these  struc¬ 
turally  more  perfect  samples.*  Additional  ^R  measure¬ 
ments  on  these  samples  are  presently  under  way. 

3.  Longitudinal  relaxation  rate,  A—l/T, 

The  temperature  dependence  of  the  relaxation  rate  of  the 
nonoscillating  component  of  the  ZF-/rSR  signal,  A(r) 
~l/r,  is  shown  in  Fig.  8(c).  As  a  function  of  decreasing 
temperature  A.{T)  can  be  seen  to  exhibit  a  cusplike  feature  at 
the  ferromagnetic  transition  of  the  Ru  moments  at  Tc=  133 
K  and  a  steplike  increase  at  very  low  temperature  which 
most  likely  is  related  to  the  ordering  of  the  Gd  moments.  The 
cusp  feature  at  Tc=  133  K  characterizes  the  slowing  down  of 
the  spin  dynamics  of  the  Ru  moments  as  the  ferromagnetic 
transition  is  approached.  The  cusp  maximum  occurs  when 
the  spin-fluctuation  rate  equals  &e  typical  /rSR  time  scale 
for  T^-IO"*."  Note,  that  in  the  ferromagnetically  ordered 
state  that  longitudinal  relaxation  rate  remains  unusually  large 
with  values  of  A(T<tTc)^0.2-0A  /is"^  that  are  at  least  an 
order  of  magnitude  larger  than  expected  for  a  classical 
ferromagne^^  (where  two-magnon  excitations  provide  the 
major  contribution  to  spin  dynamics).  We  have  confirmed  by 
a  fJtSR  measurements  in  a  longitudinal  field  of  H^=  6  kOe 
that  this  large  relaxation  rate  is  indeed  characteristic  for  the 
longitudinal  component  of  the  fiSR  signal.  At  present  we 
cannot  provide  a  definite  explanation  of  the  origin  of  the 
unusually  large  value  of  A.  However,  we  empha^e  that  the 
RuSr2GdCu208  system  can  be  expected  to  exhibit  a  rathw 
complex  magnetic  behavior  since,  besides  the  ferromagneti¬ 
cally  ordered  Ru  moments,  it  also  contains  the  larger  Gd 
moments  with  /i(Gd^‘")«7.4/tB  which  remain  paramagnetic 
below  Tc.  The  magnetic  ordering  transition  of  the  Gd  mo¬ 
ments  at  7’«2.6  K  is  evident  in  the  ZF-/iSR  data  in  Figs. 
8(a)-8(c)  firom  the  sudden  increase  in  the  local  magnetic 
field  (or  tile  ^ttSR  precession  firequency,  (v^))  and  a  corre¬ 
sponding  increase  in  botii  relaxation  rates,  k  and  A.  In  addi¬ 
tion,  we  note  tiiat  recently  it  has  been  shown  by  /tSR  mea¬ 
surements  that  in  strongly  underdoped  high-TcCuprate 
superconductors  (like  the  present  RuSr2GdCu208  compound) 
also  the  Cu  moments  exhibit  a  spin-glass-type  fisezing  traii- 
sition  at  low  temperature.^  Finally,  it  appears  that  the  longi¬ 
tudinal  relaxation  rate  A  exhibits  m  ^ditional  weak 
anomaly  at  a  temperature  of  r**20  K,  i.e.,  in  the  vicinity  of 
the  superconducting  transition  at  Tc  —  16  K.  At  present  we 
are  not  sure  whether  this  effect  is  related  to  the  onset  of 
superconductivity.  From  Fig.  8(b)  it  appears  that  the  ttans- 
verse  relaxation  rate  also  exhibits  a  steplike  increase  in  the 
same  temperature  range.  The  local  magnetic  field  at  the 
muon  site,  however,  [see  Fig.  8(a)]  does  not  seem  to 
exhibit  any  anomaly  in  the  vicinity  of  Tg.  We  expect  that 
further  /rSR  measurements  on  rare-earth  substituted 
RuSr2Gd,_;ci?*Cu20*  samples,  as  weU  as  on  less  strongly 
underdoped  samples  with  higher  critical  temperatures  of 


TABLE  I.  The  local  magnetic  field  at  the  miion  site  (B^),  ob¬ 
tained  fi^om  the  dipolar-field  calculation.  The  results  are  shown  for 
two  different  muon  sites  and  for  three  different  orientations  with  the 
Ru  moments  [/i(Ru*^)=  I /Xj,]  ferromagnetically  ordered  along  the 
Ru-0  bond  [100],  along  the  diagonal  [1 10],  or  perpendicular  to  the 
Ru02  planes  [001],  Note  that  for  the  [110]  orientation  there  exist 
two  magnetically  inequivalent  muon  sites. 


Muon  site 

[100] 

[110] 

[001] 

(0.13,0.225,0.16) 

(0.128,0.222,0.175) 

960  G 
740  G 

1068/805  G 
824/665  G 

1471  G 
1231  G 

up  to  40  K,  should  shed  more  light  on  the  complex  magnetic 
behavior  and  its  interplay  with  superconductivity  in  the  Ru- 
1212  system. 

C.  Dipolar  field  calculation 

While  the  ZF-/xSR  data  give  clear  evidence  for  the  pres¬ 
ence  of  a  homogeneous  magnetically  ordered  state,  they  do 
not  provide  any  direct  information  about  the  origin  of  the 
magnetic  moments,  the  type  of  the  magnetic  order,  and  its 
direction.  Based  on  dipolar-field  calculations  of  the  local 
magnetic  field  at  the  muon  site,  however,  one  can  test  the 
consistency  with  an  assumed  magnetic  structure.  The  result 
of  these  calculations  depends  on  the  location  of  the  intersti¬ 
tial  muon  site  and  also  on  the  orientation  of  the  Ru  moments. 
Unfortunately,  for  the  Ru-1212  system  neither  of  these  is 
accurately  known  at  present.  Nevertheless,  it  seems  plausible 
that  the  muon  site  is  similar  to  that  in  YBa2Cu307_f  (and 
other  related  cuprate  compounds)  where  the  positive  muon 
forms  a  hydroxyl  bond  with  the  apex  oxygen  and  is  located 
at  the  so-called  “apical  site”  near  the  point 
(0.12a,6.225Z>,0.14c).‘^  Indeed,  as  is  summarized  in  Table  I, 
we  obtain  rather  good  agreement  with  the  experimental  value 
of  (S  )(r->0)  =  720  G  if  we  take  a  similar  apical  site  near 
the  point  (0.13a,  0.226,  0.1 6-0. 17c)  and  assume  that  the 
ferromagnetically  ordered  Ru  moments  [/i(Ru)=l/rB]^re 
oriented  along  tiie  Ru02  plane  either  along  the  Ru-0  bond, 
[100],  or  along  the  diagonal  [110]  (see  Table  I).  For  the 
[110]  orientation,  however,  there  exist  two  magnetically  in¬ 
equivalent  muon  sites  which  should  give  rise  to  two  distinc¬ 
tive  precession  frequencies  in  the  ;xSR  spectra  (which  are  not 
observed  experimentally).  For  the  Ru  moments  oriented  per¬ 
pendicular  to  the  Ru02  layer  along  [001]  the  resulting  local 
magnetic  field  at  the  apex  site  is  significantly  larger  than  the 
experimental  value.  In  order  to  obtain  reasonable  agreement 
with  experiment  for  the  [001]  orientation,  one  has  to  assume 
that  the  muon  site  is  located  much  closer  to  the  Cu02  planes. 
Such  a  muon  site,  however,  is  not  very  realistic  (^ply 
speaking  the  positive  muon  is  repelled  by  the  positively 
charged  Cu02  planes)  and  has  not  been  obse^ed  in  any  of 
the  related  cuprate  compounds.  We  thus  tentatively  conclude 
that  the  moments  align  in-plane  consistent  with  the  two- 
component  XY  scenario.  While  this  result  is  rather  conve¬ 
nient  in  terms  of  the  coexistence  of  the  ferromagnetic  order 
of  the  Ru  moments  and  the  superconductive  which  resides 
within  the  Cu02  layers  as  discussed  below,  one  has  to  keep 
in  that  the  underlying  assumptions  are  rather  crude.  For 

more  detailed  and  decisive  information  on  the  structure  and 
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the  orientation  of  the  Ru  spin  order  we  must  await  the  result 
of  neutron-scattering  experiments. 

IV.  A  POSSIBLE  SCENARIO  FOR  COEXISTENCE  OF 
FERROMAGNETIC  AND  SUPERCONDUCTING  ORDER 

Having  established  that  the  ferromagnetic  and  the  super¬ 
conducting  order  parameter  coexist  on  a  microscopic  scale, 
we  arrive  at  the  important  question  as  to  how  this  system 
manages  to  avoid  strong  pair-breaking  effects.  We  suspect 
that  the  answer  is  closely  related  to  the  layered  structure  of 
the  hybrid  ruthenate-cuprate  compound  and,  in  particular,  to 
the  purely  two-dimensional  coherent  charge  transport  in  the 
strongly  underdoped  Cu02  planes.  We  envisage  a  scenario 
where  the  ferromagnetically  ordered  Ru  spins  are  aligned  in 
the  Ru02  plane  having  a  very  large  out-of-plane  anisotropy 
while  the  charge  dynamics  of  the  superconducting  Cu02 
planes  is  purely  two-dimensional,  i.e.,  coherent  charge  trans¬ 
port  occurs  only  along  the  direction  of  the  Cu02  planes.  For 
such  a  configuration  the  principal  pair-breaking  effect  due  to 
the  electromagnetic  interaction  can  be  minimized,  since  the 
dot  product  of  the  magnetic  vector  potential  (which  then  is 
normal  to  the  planes)  and  the  momentum  of  the  Cooper  pair 
(which  is  parallel  to  the  planes)  vanishes.  An  additional  re¬ 
quirement  is  that  the  direct  hyperfine  interaction  between  the 
superconducting  electrons  of  the  Cu02  planes  and  the  or¬ 
dered  Ru  spins  has  to  be  extremely  small.  Both  requirements 
may  be  fulfilled  in  the  present  Ru-1212  system  due  to  the 
confinement  pf  the  superconducting  electrons  of  the  steongly 
underdoped  Cu02  planes.  The  absence  of  magnetic  pair 
breaking  is  suggested  by  the  fact  that  the  value  is  fully 
consistent  with  the  underdoped  state  indicated  by  the  thp- 
moelectric  power.  So  far  we  have  been  unable  to  rignifi- 
cantly  increase  the  doping  state  (by,  e.g.,  Ca  substitution)  so 
as  to  explore  these  implications.  Furthermore,  we  have  not 
yet  succeeded  in  crystallographically  aligning  powdem  or 
growing  single  crystals  which  would  allow  one  to  investigate 
the  magnetic  anisotropy.  However,  further  information  can 
be  obtained  by  examining  the  ferromagnetic  fluctuations 
above  Tq  as  seen  in  the  divergence  of  the  susceptibility. 
Figure  9  shows  d^ldT  plotted  versus  (7Vrc““l)  for  the 
zero-field-cooled  susceptibility  for  T>  T q  .  The  slope  of 
-2.30(3)  indicates  a  critical  exponent  of  1.30(3)  consis¬ 
tent  with  3D  JllT  fluctuations  for  which  y=1.32  (Ref.  18)  and 
again  consistent  with  orientation  of  the  Ru  moments  within 
the  a-b  plane. 

Finally,  we  note  that  an  alternative  (and  highly  specula¬ 
tive)  explanation  for  the  coexistence  of  high-T^.  supercon¬ 
ductivity  and  ferromagnetic  order  in  the  present  Ru  1212 
superconductor  could  be  that  the  superconducting  order  pa¬ 
rameter  has  a  nonzero  angular  momentum  which  itself 
breaks  time-reversal  symmetry.  Such  an  unconventional  or¬ 
der  parameter  symmetry  has  been  discussed  also  in  the  con¬ 
text  of  the  Sr2Ru04  superconductor.  We  point  out,  however, 
that  at  present  we  have  no  evidence  in  favor  of  such  a  sce¬ 
nario. 


FIG.  9.  d^/dT  plotted  versus  {T/Tc-l)  for  the  zero-field- 
cooled  susceptibility  in  the  temperature  range  of  T>  Tc-  The  slope 
of  -2.30(3)  indicates  a  critical  exponent  of  y=  1.30(3)  consistent 
with  3D  A'r  fluctuations  for  which  y=1.32  (Ref.  18)  and  consistent 
with  orientation  of  the  Ru  moments  within  the  a^b  plane. 

V.  CONCLUSIONS 

In  summary,  we  have  performed  dc-magnetization  and 
zero-field  muon-spin  rotation  (ZF-/iSR)  measurement 
which  characterize  the  superconducting  and  the  magnetic 
properties  of  the  hybrid  cuprate-ruthenate  compound 
RuSr2GdCu208.  The  dc-magnetization  data  establish  that  this 
material  exhibits  ferromagnetic  order  (or  at  least  magnetic 
order  with  a  sizeable  ferromagnetic  component)  below  Tc 
=  133  K  and  becomes  superconducting  at  a  much  lower  tem¬ 
perature  of  r^=  1 6  K.  We  obtain  evidence  that  superconduct¬ 
ing  charge  carriers  originate  firom  the  Cu02  planes,  while  the 
ferromagnetic  order  is  associated  with  the  Ru  moments  with 
/x(Ru)«  1  fJiB .  The  larger  Gd  moments  with  /i(Gd)«=‘7.4^5 
do  not  appear  to  participate  in  the  ferromagnetic  transition 
but  remain  paramagnetic  to  very  low  temperature  and  un¬ 
dergo  most  likely  an  antiferromagnetic  transition  at 
=  2.6  K.  The  ZF-/xSR  experiments  provide  evidence  that  the 
ferromagnetic  phase  is  homogeneous  on  a  microscope  scale 
and  accounts  for  most  of  the  sample  volume.  Furthermore, 
they  indicate  that  the  magnetically  ordered  state  is  not  sig¬ 
nificantly  modified  by  the  onset  of  superconductivity.  This 
rather  surprising  result  raises  the  question  as  to  how  ferro¬ 
magnetic  and  superconducting  order  can  coexist  on  a  micro¬ 
scopic  scale  while  avoiding  strong  pair-breaking  effects  that 
tend  to  destroy  superconductivity.  We  have  outlined  a  pos¬ 
sible  scenario  which  relies  on  the  two-dimensional  charge 
dynamics  of  the  Cu02  planes  and  the  assumption  that  the 
ferromagnetic  order  parameter  of  the  Ru  moments  is  con¬ 
fined  to  the  Ru02  layers. 
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Abstract.  Muon  spin  relaxation,  as  a  local  magnetic  probe,  is  not  normally  expected  to  be  able  to 
provide  measurements  of  any  magnetic  correlation  length.  When  anomalously  shallow  static  zero- 
field  muon  spin  relaxation  functions  were  observed  for  CeNii^-jcCu^Sn  and  a  small  number  of  other 
highly  disordered  magnetic  systems,  however,  and  a  ‘Gaussian-broadened  Gaussian*  relaxation 
function  shown  to  fit  them,  it  was  suggested  that  an  unusual  short-range  correlation  among  local 
moments  in  the  material  was  acting.  ‘Range-correlated  moment  magnitude  variation’,  in  which 
ion-moment  orientations  are  completely  random,  but  moment  magnitude  varies  slowly  from  place 
to  place,  is  shown  by  Monte  Carlo  numerical  simulation  on  a  lattice  of  moments  to  generate  the 
qualitative  form  of  shallow  static  zero-field  relaxation  observed.  Thus,  the  shape  of  the  zero-field 
muon  spin  relaxation  function  is  sensitive  to  a  correlation  length. 


1.  Introduction 

Positive-muon  spin  relaxation  (f^SR)  is  a  probe  of  local  magnetic  fields  in  materials  (for 
introductory  reviews,  see  [1])  that  is  particularly  useful  when  carried  out  in  zero  external  field 
(ZF).  The  magnetic  moment  of  a  fi*  stopped  at  an  interstitial  site  in  a  solid  processes  in  the  local 
field  at  that  site,  and  then  emits  a  positron  preferentially  in  the  direction  of  the  /^'^-moment  at  the 
instant  of  decay.  Knowing  the  initial  polarization  of  the  muons  and  detecting  the  directions  of 
the  decay  positrons,  the  experimenter  can  measure  the  time  evolution  of  the  polarization  in  the 
direction  of  initial  orientation,  called  the  (longitudinal)  ‘muon  spin  relaxation  function’,  G^Ct). 
In  ferromagnetic  or  antiferromagnetic  states,  the  well-ordered  array  of  ion  moments  will  often 
generate  a  unique  field  at  the  muon  site,  causing  a  single  muon  precession  firequency  observed 
as  spontaneous  oscillation  of  G^(t)  (for  reviews  of  fiSR  in  magnetic  materials,  see  [2]).  In 
cases  where  the  magnetic  moments  in  the  material  are  static  (not  fluctuating)  yet  not  long- 
range  ordered  (as  in  frozen  spin  glasses  and  conductors  containing  only  nuclear  moments), 
the  ‘Kubo-Toyabe’  relaxation  functions  [3]  corresponding  to  either  a  Gaussian  local  field 
distribution 

Gf(/)  =  5.5(l-AV).xp(-^)  (1) 

(where  A  =  y^(B,)rmj.  and  is  the  muon  gyromagnetic  ratio)  or  to  a  Lx)rentzian  local  field 
distribution 

G^(t)  =  ^  exp(-aO  (2) 

(where  a  =  Y^{Bx)hwhm)  have  been  sufficient  to  fit  all  data  obtained  until  recently.  In  both 
of  these  functions,  polarization  drops  from  the  initial  value  of  unity  to  a  single  mii^um,  and 
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then  recovers  to  1/3  at  late  times.  This  1/3  asymptote  is  a  characteristic  feature  of  muon  spin 
relaxation  in  static  random  local  fields  and  zero  external  field.  It  is  now  known,  however, 
that  in  a  small  number  of  materials,  including  Al-Mn-Si  magnetic  quasicrystals  [4]  and 
CeNii-jcCu;cSn  Kondo-lattice  alloys  [5],  static  ZF  relaxation  of  initial  polarization  can  lead  to 
a  minimum  too  shallow  to  be  reproduced  by  either  of  the  Kubo-Toyabe  functions  (equation  (1) 
or  equation  (2))  above.  In  fact,  the  extreme  case  of  monotonic  decay  to  the  1/3  asymptote  has 
been  observed. 

A  prior  publication  [6]  showed  that  the  observed  relaxation  in  CeNii-.;cCujcSn  alloys  is 
fitted  well  by  a  relaxation  function  derived  from  a  ‘convolution  broadening’  of  a  Gaussian 
distribution.  This  means  that  instead  of  a  single-Gaussian  random  local  field  distribution  with 
a  unique  value  of  {Bx)rms^  a  collection  of  distributions  are  used,  corresponding  to  a  collection 
of  different  muon  sites,  where  the  collection  of  values  of  {Bx)rms  themselves  form  a  Gaussian 
distribution  of  width  W  and  most  likely  value  B^.  This  ‘Gaussian-broadened  Gaussian’  static 
ZF  Kubo-Toyabe  relaxation  function  has  the  form 


qGBG 


(0  = 


1^2/ 

3  3V1  +  /?2  +  /?2a2^/V  \  \  +  + 


X  expl 


,2(1  + 


(3) 


where  and  R  =  W/Bq.  Previously,  however,  no  microscopic  model  of 

moments  in  a  solid  was  presented  that  could  produce  such  a  field  distribution. 

Directly,  /xSR  ‘senses’  only  the  magnetic  fields  at  the  muon  sites,  and  as  such  is  a  local 
probe  that  has  not  been  expected  to  be  sensitive  to  correlation  lengths.  In  cases  that  require 
a  non»trivial  value  of  R  in  equation  (3),  however,  that  non-trivial  value  (or,  equivalently, 
the  shallowness  of  the  minimum  of  polarization)  is  extra  information  not  provided  by  the 
standard  Gaussian  or  Lorentzian  forms.  One  publication  [4]  had  speculated  that  this  kind  of 
anomalous  relaxation  could  be  produced  by  random  variation  of  the  magnitude  of  the  frozen 
moments,  with  a  correlation  range.  The  present  paper  presents  the  results  of  Monte  Carlo 
simulations  of  the  ZF  muon  spin  relaxation  functions  generated  by  a  microscopic  model  of 
such  ‘range-correlated  moment  magnitude  variation’  (RCMMV).  This  involves  similar  values 
of  the  magnitude  of  nearby  magnetic  moments  that  nonetheless  have  uncorrelated  orientations. 
The  results  demonstrate  that  this  type  of  correlation  in  a  magnetically  frozen  (but  disordered) 
state  will  generate  the  kind  of  shallow  static  ZF  Kubo-Toyabe  relaxation  function  that  has 
been  observed  for  the  materials  listed  above.  As  the  correlation  length  gets  successively 
larger,  the  associated  relaxation  function  more  closely  approaches  monotonic  relaxation  to  the 
1/3  asymptote,  as  will  be  discussed.  The  extra  information  represented  by  the  value  of  R  is 
indeed  the  range  of  this  correlation.  In  the  model  discussed  here,  this  correlation  is  imposed, 
not  generated  by  a  plausible  physical  mechanism.  The  identification  of  candidate  mechanisms 
is  a  subject  for  future  research. 

Throughout  this  discussion,  it  is  necessary  to  understand  that  while  the  Lorentzian  Kubo- 
Toyabe  relaxation  function  has  a  shallower  minimum  than  the  Gaussian  Kubo-Toyabe  function, 
G^(t)  does  not  fit  any  data  of  the  type  exemplified  by  the  CeNii-jcCu^Sn  data.  The  Lorentzian 
case  couples  that  shallower  minimum  with  rapid  initial  relaxation.  If  forced  to  fit  the  initial 
relaxation  observed  in  CeNii^xCu^Sn  [5],  the  Lorentzian  minimum  occurs  far  later  than  in  the 
data.  If  forced  to  fit  the  observed  minimum,  the  initial  decay  of  G^(t)  is  much  faster  than  that 
of  the  data.  The  Gaussian-broadened  Gaussian  function  places  the  minimum  at  essentially  the 
same  time,  relative  to  the  initial  relaxation,  as  the  standard  Gaussian  Kubo-Toyabe  function, 
and  does  fit  the  data  well.  Since  Gaussian-like  behaviour  is  generally  associated  with  dense- 
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moment  systems  [2],  a  dense  array  of  moments  was  specifically  incorporated  into  the  Monte 
Carlo  model. 

2.  The  microscopic  Monte  Carlo  model 

A  Monte  Carlo  simulation  of  a  microscopic  model  of  moments  on  a  lattice  has  been  used  before 
to  simulate  the  field  distribution  at  an  interstitial  muon  site  and  the  corresponding  static  ZF 
muon  spin  relaxation  function  [7].  For  a  cluster  of  a  finite  number  of  ions  on  a  lattice,  a  vector 
magnetic  moment  is  assigned  to  each  magnetic  ion,  the  assumed  coupling  of  the  moments  to  a 
muon  at  a  particular  site  is  used  to  calculate  the  effective  field  due  to  each  ion  at  the  muon  site, 
and  those  individual-ion  fields  are  vector  summed  to  find  the  net  local  field  at  the  muon  for  that 
particular  moment  configuration.  For  any  particular  initial  muon  polarization,  the  precession 
of  the  muon  spin  around  the  local  field  is  calculated.  Of  particular  interest  is  the  projection 
of  the  polarization  along  its  initial  direction,  as  a  function  of  time.  This  is  the  longitudinal 
relaxation  function  G^it).  For  a  single  muon,  this  will  oscillate  at  the  Larmor  frequency  of  the 
local  field.  If  there  is  any  disorder  in  the  system,  randomization  characteristic  of  that  disorder 
is  included  in  the  calculation,  and  the  calculation  is  repeated  many  times.  The  many  net  local 
fields  are  collected  into  histograms  representative  of  the  local  field  probability  distribution,  and 
the  projected  polarization  of  many  muons  is  averaged  to  generate  the  simulated  static  muon 
spin  relaxation  function. 

One  type  of  disorder  that  often  occurs  in  samples  studied  with  fiSR  is  polycrystallinity. 
Most  /zSR  experiments  are  performed  on  powders  or  polycrystalline  ingots,  where  all  jjossible 
orientations  of  the  principal  axes  of  the  crystallites  of  the  material,  with  respect  to  the  initial 
muon  polarization,  occur.  This  averages  over  any  anisotropy  of  the  relaxation  with  respect  to  the 
crystal  axes.  All  of  the  equations  above  assume  ‘polycrystalline  averaging’  (or  perfect  isotropy 
of  the  local  field  distribution).  In  a  Monte  Carlo  program,  this  is  achieved  by  randomizing  the 
orientation  of  the  lattice  principal  axes  relative  to  the  initial  muon  polarization,  and  this  was 
done  in  every  case  in  the  simulations  discussed  below. 

A  previous  publication  [7]  discussed  Monte  Carlo  simulation  of  the  static  ZF  relaxation 
functions  that  should  be  generated  by  non-dilute  alloys,  and  that  the  path  to  the  Lorentzian 
case  in  the  dilute  limit  was  not  simple.  That  program  used  fixed-magnitude,  random-direction 
moments  in  a  cluster  of  lattice  positions.  Upon  discovery  of  the  anomalous,  shallow-minima, 
static  ZF  relaxation  functions,  and  with  the  expectation  that  some  new  type  of  disorder 
would  be  required  in  the  explanation,  a  simple  modification  was  made  to  that  program  which 
randomized  the  moment  magnitude  in  addition  to  the  orientation.  The  modification  produced 
simulated  field  distributions  and  relaxation  functions  of  the  same  qualitative  form  as  those 
shown  in  the  earlier  paper.  They  evolved  from  the  Gaussian  Kubo-Toyabe  function  at  high 
moment  concentration,  through  two-minima  forms  (with  two  effective  sites)  at  intermediate 
concentrations,  to  the  Lorentzian  limit  at  low  concentrations,  without  producing  anything  that 
fitted  the  shallow  relaxation  functions  of  CeNii_;tCu;tSn  and  magnetic  Al-Mn-Si  quasicrystals. 
Completely  random  moment  variation  does  not  explain  the  anomalous  relaxation  that  has  been 
observed. 

The  subject  of  this  paper  is  a  model  that  more  successfully  reproduces  shallow  static  ZF 
relaxation  functions.  It  was  constructed  to  be  a  specific  example  of  range-correlated  moment 
magnitude  variation,  alluded  to  above.  Retaining  a  lattice  cluster  and  completely  random 
moment  orientations,  moment  magnitudes  are  held  fixed  in  sub-clusters  of  moments,  while 
they  are  varied  randomly  between  sub-clusters.  Forcing  the  sub-clusters  to  have  a  particular 
average  size  creates  a  correlation  range.  To  obtain  good  statistical  behaviour  of  the  simulations 
using  such  sub-clusters,  substantially  larger  overall  cluster  size  is  necessary:  while  the  clusters 
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of  the  previous  Monte  Carlo  fiSR  paper  [7]  were  spheres  of  up  to  500  ions  on  fee  lattice  sites, 
the  work  reported  here  used  cubes  of  ions  on  simple-cubic  lattice  sites,  with  up  to  126  ions 
on  a  side  (containing,  therefore,  up  to  slightly  more  than  two  million  ions).  The  use  of  cubic 
clusters,  and  a  simple-cubic  lattice,  help  reduce  CPU  time,  as  does  the  increase  in  computation 
speed  available  since  that  previous  publication.  The  qualitative  conclusions  drawn  should  not 
depend  on  the  particular  lattice  used,  although  fine  details  may  well  do  so. 

Within  an  overall  cluster,  the  sub-clusters  were  chosen  to  be  rectangular  boxes,  with  each 
box’s  three  side  lengths  chosen  randomly  between  one  ion  thick  and  a  maximum  thickness  (the 
‘cut-off’,  which  controls  the  correlation  range),  subject  to  a  simple  packing  constraint  such 
that,  in  spite  of  the  randomization,  the  overall  cube  is  filled  by  the  boxes,  with  no  ions  left 
out.  For  each  box,  a  single  value  of  the  moment  magnitude  is  chosen  randomly  between  zero 
and  unity.  This  rectangular-box  configuration  of  correlated  clusters  is  not  physically  likely, 
but  is  computationally  simple.  More  realistic  methods  of  generating  correlations,  usually 
involving  some  pairwise  nearest-neighbour  interaction,  would  require  orders  of  magnitude 
more  CPU  time,  which  would  be  prohibitive.  Furthermore,  no  physically  reasonable  candidate 
mechanisms  for  producing  such  correlations  are  known  at  this  time.  The  simulations  reported 
here  provide  a  first  glimpse  of  what  RCMMV  can  do  to  static  ZF  relaxation  functions. 

More  challenging,  in  terms  of  the  microscopic  physics,  is  the  assumption  of  continuous 
moment  variation  between  zero  and  a  maximum.  Ion  magnetic  moments  in  solids  are  usually 
single  valued,  and  cases  of  static  moment  instability  are  usually  discussed  in  terms  of  (usually 
only  two)  distinct  values  (‘mixed  valence’ ;  see,  e.g.,  [8]).  A  continuous  range,  or  at  least  a  large 
number  of  possible  values,  of  moment  magnitude  seems  necessary  to  generate  shallow  static  ZF 
relaxation.  Use  of  the  models  described  here  with  only  two  possible  moment  values  generates 
field  distributions  and  relaxation  functions  clearly  representative  of  two  muon  sites  with  two 
particular  rms  field  values,  and  does  not  reproduce  the  shallow  static  ZF  relaxation  observed. 
The  behaviour  discussed  here  might  well  be  generated  by  the  size  of  the  magnetic  ion  moment 
taking  a  distinct  and  unique  value  (not  just  on/off)  for  each  and  every  different  near-neighbour 
environment  that  occurs  in  the  material.  For  a  site-substitution  alloy  like  CeNii-.xCu;cSn,  this 
would  result  in  a  binomial  distribution  of  moment  sizes.  Continuous  ranges  of  moment  values 
are  a  possibility  in  Kondo  screening  (see,  e.g.,  [9]),  but  static  spatial  inhomogeneity  with  a 
continuous  range  of  moment  values  may  not  yet  have  been  considered  in  a  Kondo-interaction 
calculation. 

The  correlation  generated  by  the  modelling  should  be  monitored  by  a  two-point  correlation 
function,  which,  if  it  decays  exponentially  in  the  distance  r  between  the  two  points,  defines 
a  correlation  length.  Any  moment-moment  correlation  can  only  be  defined  at  those  discrete 
values  of  r  that  represent  the  separations  of  pairs  of  moments  in  the  lattice  cluster.  For  any 
such  r,  there  will  usually  be  a  number  N(r)  of  pairs  of  moments,  with  magnitudes  m,  and  nij, 
for  which  nj  =  r.  The  correlation  function  chosen  was 

•  ^*j 

At  a  particular  r,  the  sum  is  over  all  pairs  for  which  =  r.  In  this  notation, 

N(r)  =  X;(l)r„=r. 

u 

C(r)  is  unity  when  all  moments  are  the  same  size,  and  2/3  when  the  magnitudes  are  completely 
random  on  the  unit  interval.  Figure  1  shows  C(r)  for  a  typical  simulation  run,  with  an  overall 
cluster  thickness  of  44  ions  and  a  sub-cluster  thickness  cut-off  of  10  ions.  All  distances 
are  stated  as  multiples  of  the  distance  between  nearest-neighbour  moments  in  the  lattice 
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Figure  1.  The  moment  magnitude  correlation  C(r)  (random  =  2/3),  as  a  function  of  distance  r 
between  the  moments  (in  units  of  the  distance  between  nearest-neighbour  moments),  for  a  typical 
Monte  Carlo  simulation  of  range-correlated  magnetic  moment  variation  (RCMMV,  as  described 
in  the  text),  with  a  cubic  cluster  thickness  of  44  ions,  and  a  rectangular-box  sub-cluster  cut-ojff 
thickness  of  10  ions.  The  solid  curve  is  a  fit  of  exponential  decay  to  the  baseline,  with  a  deduced 
correlation  length  of  3.29(5). 


Figure  2.  The  RCMMV  moment  magnitude  correlation  length  deduced  from  fits  to  C(r)  as  a 
function  of  the  sub-cluster  cut-off  thickness.  The  straight  line  (slope  0.29)  is  a  guide  for  the  eye. 


(abbreviated  as  nn).  It  was  found  that  for  runs  with  cut-off  greater  than  ^1/4  of  the  overall 
thickness,  the  correlation  did  not  settle  down  near  2/3  for  large  r,  so  these  were  not  used.  Also 
shown  in  figure  1  is  the  fit  of  exponential  decay  to  the  baseline, 

C(r)  =  +  B  (5) 

(where  A  1/3  and  B  ~  2/3),  with  a  deduced  correlation  length  Xc  =  3.29  ±  0.05  nn. 
Figure  2  shows  the  correlation  length  deduced  from  such  fits  as  a  function  of  the  cut-off 
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thickness.  The  correlation  length  is  :^0.3  times  the  cut-off  thickness.  With  overall  clusters  of 
up  to  two  million  ions,  correlation  lengths  only  up  to  about  ten  times  the  nn  distance  can  be 
reasonably  modelled. 

For  eafch  overall  cluster,  once  moment  magnitudes  were  assigned  as  above,  and  moment 
directions  randomly,  numerous  examples  of  the  interstitial  muon  site  (at  the  crystallo- 
graphically  unique  unit-cell  centre)  were  chosen  at  random  from  the  interior  volume  at  least 
five  nn  steps  from  the  outside  of  the  cluster.  Modelling  of  the  field  distributions  and  static  ^ 
relaxation  functions  was  then  performed  as  described  previously  [7].  The  local  field  at  each 
muon  site  was  calculated  in  dipole  coupling  to  all  the  moments  in  the  cluster,  and  then  added  to 
the  field  distribution  histograms.  The  component  of  muon  polarization  in  its  initial  direction, 
as  a  function  of  time,  was  calculated  in  that  local  field  and  added  to  the  relaxation-function 
histogram.  For  the  larger  cluster  sizes,  typically  moments  were  assigned  600  times  and  more 
than  2000  muon  sites  were  evaluated  for  each  of  those  moment  configurations. 


-60  -40  -20  0  20  40  60 

Field  Cartesian  Component  (kG) 

Figure  3.  Monte  Carlo  RCMMV  local  field  Cartesian-component  distribution  histograms:  (a)  for 
an  overall  thickness  of  44  ions  and  a  correlation  cut-off  thickness  of  10  ions;  the  solid  (dashed) 
curve  shows  the  Gaussian  (Lorentzian)  distribution  with  the  same  maximum  and  fwhm;  (b)  for  an 
overall  thickness  of  126  ions  and  a  correlation  cut-off  of  30  ions. 


3.  Results 

The  Gaussian  and  Lorentzian  Kubo-Toyabe  relaxation  functions  are  named  after  the  distrib¬ 
utions  of  one  Cartesian  component  of  the  magnetic  field,  P{Bx),  that  generate  them.  As 
previously  discussed,  however  [7],  particularly  in  the  case  of  poly  crystalline  averaging,  the 
distribution  of  the  magnetic  field  magnitude  PQB])  is  more  directly  related  to  the  static  ZF 
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Figure  4.  Monte  Carlo  RCMMV  local  field  magnitude  distribution  histograms:  (a)  for  an  overall 
thickness  of  44  ions  and  a  correlation  cut-off  thickness  of  10  ions;  the  solid  (dashed)  curve  shows  the 
Maxwellian  (Lorentzian  [7, 10])  magnitude  distribution  with  same  maximum  value  and  position; 
(b)  for  an  overall  thickness  of  126  ions  and  a  correlation  cut-off  of  30  ions. 


relaxation  function,  by 

Gz(f)  =  ^  +  jjf  />(|B|)cos(y^|B|f)ci|B|.  (6) 

Therefore  both  the  component  and  magnitude  distributions  were  monitored  by  the  Monte  Carlo 
program.  Examples  are  shown  in  figure  3  and  figure  4,  respectively.  For  comparison,  curves  in 
figure  3  show  the  Gaussian  and  Lorentzian  distributions  with  the  same  fwhm  as  the  simulation: 
the  RCMMV  model  generates  distributions  more  sharply  peaked  at  low  field  than  either  of  the 
standard  ones.  The  curves  in  figure  4  show  Maxwellian  and  Lorentzian  [7, 10]  magnitude 
distributions,  and  the  model  clearly  generates  distributions  with  more  probability  of  the 
occurrence  of  low  fields,  relative  to  the  most  likely  field.  Figure  5  shows  a  sequence  of  simulated 
static  ZF  relaxation  functions,  for  the  correlation  lengths  shown.  The  simulation  shown  for 
no  correlation  (every  moment  magnitude  and  direction  random)  is  almost  indistinguishable 
from  the  standard  Gaussian  Kubo-Toyabe  function.  The  model,  as  described,  naturally  keeps 
the  average  moment  (magnitude)  per  unit  volume  constant,  independently  of  the  correlation 
length,  and  this  results  in  a  constant  rms  local  field  Cartesian  component  {Bx)rms  •  For  moments 
in  the  interval  (0, 1  Mb)  at  1 A  spacing,  as  shown  for  figures  3  and  4,  -  19.8±0.2kG. 

As  a  further  result,  it  can  be  seen  in  figure  5  that  the  initial  drop  of  polarization  firom  1.0  to 
"^0.6  is  also  essentially  independent  of  the  correlation,  and  that  the  time  at  which  minimum 
polarization  occurs  is  not  strongly  dependent  on  the  correlation,  but  that  the  value  of  minimum 
polarization  achieved  increases  with  Figure  6  shows  the  minimum  polarization  of  the 
Monte  Carlo  relaxation  functions,  plotted  as  a  function  of  1  /X^. 
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Figure  5.  Monte  Carlo  RCMMV  static  ZF  muon  spin  relaxation  functions  for  the  moment  mag¬ 
nitude  correlation  lengths  indicated. 


Figure  6.  The  minimum  polarization  achieved  by  the  Monte  Carlo  RCMMV  static  ZF  muon  spin 
relaxation  functions,  as  a  function  of  the  reciprocal  of  the  moment  magnitude  correlation  length  Xc. 
The  horizontal  line  represents  the  1/3  asymptote,  above  which  the  minimum  polarization  cannot 
rise.  The  dashed  line  is  a  guide  for  the  eye. 


'Comparison  of  the  simulated  relaxation  with  the  Gaussian-broadened  Gaussian  analytic 
form  (equation  (3))  can  be  made  by  treating  each  Monte  Carlo  spectrum  as  data  and  fitting 
to  it,  as  shown  in  figure  7.  The  fits  are  reasonable,  but  not  perfect.  Figure  8  shows 
that  the  deduced  GBG  ratio  parameter  R  appears  to  have  a  simple  straight-line  relationship 
with  the  reciprocal  of  the  correlation  length,  l/kc. 

The  RCMMV  simulated  relaxation  cannot  be  properly  least-squares  fitted  to  experimental 
data,  because  only  discrete  values  of  the  cut-off  (and  hence,  the  correlation  length)  are  available. 
Figure  9  shows  a  comparison  of  simulated  RCMMV  for  Xc  =  2.94  nn  and  low-temperature 
ZF  fiSR  in  CeCuo.2Nio.8Sn  [5,6].  The  agreement  between  the  two  is  fairly  good. 
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Time  (arb) 

Figure  7.  The  Monte  Carlo  RCMMV  static  ZF  muon  spin  relaxation  function  for  a  cluster 
overall  thickness  of  44  ions  and  a  correlation  cut-off  thickness  of  iO  ions  (circles),  with  a  fit 
of  a  Gaussian-broadened  Gaussian  relaxation  function  (curve).  The  deduced  GBG  ratio  parameter 
is  /?  =  0.543(6). 


1/[Correlation  length  (nn  units)] 

Figure  8.  The  GBG  ratio  parameter  /?,  deduced  from  fits  to  Monte  Carlo  RCMMV  static  ZF  fiSR 
spectra,  as  a  function  of  1  /Xc . 


4.  Discussion 
> 

The  Monte  Carlo  model  is  crude  in  several  respects,  one  of  which  is  the  forced  total  lack 
of  correlation  beyond  the  sub-cluster  thickness  cut-off  parameter  value.  Thus  (figure  1)  it 
can  only  approximate  the  exponential  decay  of  correlation  with  distance  to  be  expected  in  a 
real  situation.  Nonetheless  it  shows  that  RCMMV  can  generate  arbitrarily  shallow  static  ZF 
relaxation  functions.  There  is  at  this  time  no  other  microscopic  model  known  to  do  so. 

The  field  distributions  generated  by  the  RCMMV  Monte  Carlo  model  have  dramatically 
enhanced  probabilities  of  low  fields  relative  to  those  of  the  average  (or  rms)  fields  of  the 
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Figure  9.  Comparison  of  the  Monte  Carlo  RCMMV  relaxation  function  having  correlation  length 
Xc  =  2.94  nn  to  the  low-temperature  ZF  /iSR  observed  in  CeCuo.aNio.gSn. 


Standard  field  distributions.  As  illustrated  in  figure  3,  P(Bx)  becomes  progressively  more 
sharply  peaked  around  5^  =  0  as  the  correlation  length  increases.  In  figure  4,  meanwhile,  the 
most  likely  field  magnitude  shifts  to  progressively  lower  values  relative  to  the  high-field  tail  of 
P(|B1).  A  strictly  phenomenological  form  which  roughly  fits  the  magnitude  distributions  is 

where  W  is  a  width  parameter.  The  power  p  is  close  to  2  for  Xc  =  0  (where  the  simulation 
magnitude  distribution  is  almost  indistinguishable  for  the  Maxwellian  case,  and  W  =  {Bx)rms)y 
but  decreases  as  Xc  increases.  A  graph  of  p  as  a  function  of  1  /Xc  suggests  that  in  the  limit  of 
long  correlation  length,  p  0;  that  is,  the  magnitude  distribution  would  be  approximately 
Gaussian  (with  W  =  IBlrms)-  No  such  closed-form  approximation  has  yet  been  found  for 
P(Bx). 

Note  that  the  relaxation  function  Gz(t),  as  determined  from  equation  (3),  is  (a  constant 
plus)  the  cosine  Fourier  transform  of  P{\B\).  From  this,  the  (extrapolated)  change  of  the 
magnitude  distribution  from  possessing  a  peak  at  non-zero  field  to  possessing  a  peak  only 
at  |B|  =  0,  as  Xc  oo,  implies  a  change  in  from  underdamped  (possessing  a  local 
minimum)  to  overdamped  (showing  monotonic  relaxation  to  the  asymptote)  in  that  limit. 
Similarly,  the  extrapolation  of  the  minimum  polarization  on  figure  6  suggests  that  the  local 
minimum  disappears  into  the  1/3  asymptote  for  infinite  Xc.  It  thus  appears  that  the  complete 
range  of  shallow  minima,  down  to  vanishingly  small,  is  spanned  by  this  model. 

The  extrapolation  of  the  Gaussian-broadened  Gaussian  ratio  R,  in  figure  8,  however, 
suggests  a  limit  value  for  R  of  slightly  less  than  0.8,  which  still  implies  a  small  local  minimum 
of  "polarization.  The  fits  of  to  the  Monte  Carlo  relaxation  functions  are  not  perfect, 
and,  to  fit  the  Monte  Carlo  initial  relaxation  better,  overestimate  the  depth  of  the  minimum  for 
the  larger-Xc  cases.  This  discrepancy  may  be  due  to  the  crudeness  of  the  Monte  Carlo  model 
(but  that  idea  can  only  be  tested  by  developing  models  that  are  more  realistic).  It  would  be 
nice  to  be  able  to  use  figure  8  as  a  calibration  curve  for  conversion  of  measured  GBG  ratio 
values  to  underlying  correlation  lengths,  but  the  imperfect  match  of  Gf^^(t)  to  the  RCMMV 
Monte  Carlo  relaxation  functions  limits  confidence  in  this.  Direct  comparison  of  individual 
simulations  to  observed  data  is  useful,  but  is  not  as  impartial  as  least-squares  fitting  is:  there 
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is  a  risk  of  systematic  personal  bias  in  the  deduced  parameter  values.  For  the  data  of  figure  9, 
since  the  simulation  for  =  3.29  nn  looks  almost  as  good  as  the  2.94  nn  case  shown,  it  is 
reasonable  to  conclude  that  =  3.1  ±  0.3  nn. 

Care  must  be  exercised  in  thinking  about  the  statistics  of  the  model  in  the  large-correlation- 
length  limit,  ‘Xc  cx)\  In  this  limit,  C(r)  becomes  unity,  and  it  would  seem  that  all 
moments  are  the  same  size,  which  is  a  case  well  known  to  generate  a  relaxation  function 
well  approximated  by  the  Gaussian  Kubo-Toyabe  function,  not  the  monotonic  relaxation  to 
1/3  asserted  above.  The  difference  is  that  for  any  finite  r,  no  matter  how  large,  as  many  muons 
will  find  sites  surrounded  by  very  small  moments  as  will  land  at  sites  with  moments  near  the 
maximum  value.  For  all  finite  r ,  the  moments  are  uniformly  spread  over  the  unit  interval.  This 
is  only  a  formal  problem.  Physically,  a  correlation  length  Xc  30  nn  units  (typically  less  than 
200  A  in  real  materials)  would  generate  a  static  relaxation  function  with  a  minimum  so  shallow 
that  it  would  be  well  approximated  by  monotonic  relaxation,  within  the  statistical  accuracy  of 
a  normal  fiSR  spectrum,  and  thus  would  appear  to  be  at  the  limit,  yet  with  sample  dimensions 
of  at  least  tenths  of  millimetres,  there  will  still  be  room  for  millions  of  such  correlated  regions 
with  vastly  different  moment  magnitudes. 


5.  Conclusions 

The  observed  anomalous  form  (with  the  characteristic  shallow  minimum)  of  static  ZF 
muon  spin  relaxation  in  disordered  magnetic  materials  such  as  CeNii_;cCu;cSn,  which 
was  previously  found  to  conform  to  the  phenomenological  ‘Gaussian-broadened  Gaussian’ 
relaxation  function,  has  been  shown  here,  by  Monte  Carlo  simulation,  to  be  associated  with 
a  subtle  short-range  partial-ordering  effect:  range-correlated  moment  magnitude  variation 
(RCMMV).  The  shallowness  of  the  ZF  relaxation,  relative  to  standard  Gaussian  Kubo-Toyabe 
relaxation,  is  in  a  one-to-one  relationship  with  the  moment  magnitude  correlation  length,  and 
so  fiSR,  as  a  local  probe,  in  these  circumstances  provides  a  measure  of  a  correlation  length.  It 
is  a  strange  and  heretofore  unexpected  correlation,  however,  of  static  spatial  inhomogeneity  of 
the  moment  magnitudes:  the  moments  not  just  ‘on’  or  ‘off’,  but  ranging  over  enough  distinct 
values  to  be  approximated  by  a  continuous  range  from  zero  to  a  maximum. 
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